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ABSTRACT 
We analyse the distribution of CO brightness temperature and integrated intensity in M5 1 at ~ 40 pc resolution 
using new 12 CO(7 = 1 — > 0) data from the Plateau de Bure Arcsecond Whirlpool Survey (PAWS). We present 
probability distribution functions (PDFs) of the CO emission within the PAWS field-of-view, which covers 
the inner ~ 1 1 x 7 kpc of M51. We find clear variations in the shape of CO PDFs both within different M51 
environments, defined according to dynamical criteria, and between M51 and two nearby low -mass galaxies, 
M33 and the Large Magellanic Cloud. Globally, the PDFs for the inner disk of M51 can be represented by 
narrow lognormal functions that cover ~ 1 to 2 orders of magnitude in CO brightness and integrated intensity. 
The PDFs for M33 and the LMC are narrower and peak at lower CO intensities, consistent with their lower gas 
surface densities. However, the CO PDFs for different dynamical environments within the PAWS field depart 
significantly from the shape of the global distribution. The PDFs for the interarm region are approximately log- 
normal, but in the spiral arms and central region of M51, they exhibit diverse shapes with a significant excess 
of bright CO emission. The observed environmental dependence on the shape of the CO PDFs is qualitatively 
consistent with changes that would be expected if molecular gas in the spiral arms is characterised by a larger 
range of average densities, gas temperatures and velocity fluctuations, though further work is required to dis- 
entangle the relative importance of large-scale dynamical effects versus star formation feedback in regulating 
these properties. We show that the shape of the CO PDFs for different M5 1 environments is only weakly related 
to global properties of the CO emission, e.g. the total CO luminosity, but is strongly correlated with properties 
of the local giant molecular cloud (GMC) and young stellar cluster populations, including the shape of their 
mass distributions. For galaxies with strong spiral structure such as M5 1 , our results indicate that galactic-scale 
dynamical processes play a significant role in the formation and evolution of GMCs and stellar clusters. 

Subject headings: galaxies: individual (M5 1, M33, Large Magellanic Cloud) - galaxies: ISM - ISM: molecules 
- ISM: structure 
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1. INTRODUCTION 

Although the interstellar medium (ISM) represents a minor 
fraction of the baryonic matter in galaxies, it plays an im- 
portant role in their evolution, providing the raw fuel for star 
formation, receiving and then redistributing heavy elements 
created in stellar interiors, and mediating the exchange of 
matter and energy between galaxies and the intergalactic 
medium (IGM). Conditions in the ISM are influenced by 
multiple physical processes occurring across a range of 
temporal and spatial scales, including accretion of primordial 
IGM material, protostellar jets and outflows, spiral shocks, 
supernovae, and thermal instability in the diffuse atomic gas. 
As a consequence of this diverse physics, interstellar gas 
occurs in multiple phases, with temperatures, densities and 
spatial structures that span five orders of magnitude or more. 

Due to its complex hierarchical structure, quantitative 
analysis of ISM properties and dynamics is challenging. 
The intensity of emission from ISM structu res can be 
characterised in terms of a fractal index (e.g. Elmegreen & 
Falgarone 1996) or power spectrum (e.g. |Block et al.|2 010), 



while properties of the velocity field can be summarised 
using a po wer spectrum , or by the use of structure functions 
(e.g. [Brunt et"aL 2003b. The density (and column density) 
structure of the ISM is most commonly represented using a 



probability density function (PDF) (e.g. Kainulainen et al. 
|2009) . Density PDFs of the ISM are supposed to follow a 
lognormal (LN) distribution. One of the most widespread 



explanations for a LN distribution is that the expected density 
PDF result ing from a turb ulent velocity field is lognormal 
(e.g. |Padoan et al. 1997|>. Simulations of supersonically 



turbulent isothermal gas find that the width of the density 
PDF increases with the root-mean-square (rms) Mach number 
(e.g. Padoan & Nordlund 2002), and that the precise form 
of the relationship depends on the relative importance of 
compressible an d solenoidal modes in the turbulence forcing 
(Federrath et al. 2008 2010). Comparison with observations 
of the Taurus and IC5 146 molecular clouds indicates that 
turbulent driving in the ISM must contain a significant 
compressive component, or that the width of the density PDF 
depends on additional physics that is not included in the 



simulations ( Price et al.|2011) . However these models apply 
only on the scale of individual molecular clouds - where 
the assumption of isothermality is reasonable - not to whole 
galaxies. The properties of ISM turbulence should depend 
on the energy injection scale and the height of the galactic 
disk, both of which will become relevant for galactic-scale 
systems. As noted by previous authors (e.g. Wada & Norman 
2007 ), the PDFs of galaxies and individual molecular clouds 
cannot be directly compared since a galaxy's total molecular 
gas content cannot be characterised by a single temperature 
or a spatially uniform, time-independent Mach number. 

In spite of these considerations, LN PDFs are still apparent in 
galaxy-scale simulations (e.g. Wada & Norman 2007; Dobbs 
et al. 2008 Tasker & Bryan 2006J! The lognormal shape 
of the density PDFs emerges quickly in the simulations (i.e. 



within a local dynamical time), and is suprisingly robust to a 
diverse range of additional input physics (e.g. magnetic fields 
and energy feedback from stellar win ds and supernovae, 
|Wada & Norman|2001[|Dobbs et al.|20T"fl ). The characteristic 
density of high- and low-mass models is roughly invariant 
((#) ~ 1 cm" 3 ) for disks with initial mass densities spanning 
an order of magnitude (e.g. Wada & Norman 2007). The 
proposed explanation for why LN density PDFs appear to be 
generic is that galactic disks are globally stable, with a hierar- 
chical density structure that results from the action of a large 
number of stochastic, independent and nonlinear processes; 
in this case, the density PDF should evol ve towards a LN 
shape by the central limit theorem (e.g. Vazquez-Semadeni 
1994) . Notably, the only simulation in Dobbs e t al.| ( |20lf) 



that shows a strong departure from a LN PDF is the one with 
very low star formation efficiency and hence a low level of 
thermal and kinetic energy feedback into the ISM. In this 
case, the model galaxy disk does not achieve an equilibrium 
state, since most of the interstellar gas becomes quickly 
confined to dense, gravitationally bound clumps (C. Dobbs, 
priv. comm). 

More recently, LN density PDFs have gained more signifi- 
cance, as they have been used as the ba sis for explaining the 
Kenn icutt-Schmidt (KS) relation (Sch midt|[J959[ Kennicutt 



1998 1. The KS relation relates the surface density of gas to the 



star formation rate surface density according to T,sfr =A~E" s . 
On glo bal scales (i.e. averaged over entire star-forming 

(2.5 ± 0.7) x 10" 4 
100 normal 



disks) 
and n 



Kennicutt 



(1998) obtained A 



1.4 ±0.15 for a composite sample of 
and starburst galaxies. The KS relation can also be written 
using the molecular gas surface density only, in which case 
n is nearly linear (e.g. Wong & Blitz 2002; Bigie fet al.| 
[2008] [20TT) . That a LN density PDF forTEeinterstellar gas 
might natur ally yi eld the KS relation was first suggested by 
Elmegreen (2002 1, and later developed in m ore detail by 



Kravtso v (2003), Krumholz & McKee| p005) and |Wada &| 
Norman (2007). The PDF is used to find the fraction of gas 
above a given density threshold, which is assumed to collapse 
and form stars. Summing up the mass of gas at each density 
divided by the free fall time provides the star formation 
rate. This approach also forms the basi s of the universal st ar 
formation law recently proposed by |Krumholz et al!] ( |2012) . 

In spite of the potential importance of the ISM's hierarchical 
structure to the interpretation of extragalactic star formation 
laws, observational evidence to support a unive rsal LN 
density PDF for ga lactic disks remains scarce. |Gaustad| 
|& van Buren| ( |1993) used the Infrared Sky Survey Atlas to 
measure gas densities near ~ 1800 OB stars within ~ 400 pc 
of the Sun, obtaining a density PDF that resembles a LN 
function for densities in the range 0.1 < n < 10 cm" 3 . Wada 



et al. ( 2000) showed that the PDF of Hi column density in 
the Large Magellanic Cloud is approximately lognormal (at 
~ 15 pc resolution) over ~ 2 orders of magnitude. More 
recently, Berkhuijse n & Fletcher| ( |2008) derived average 
volume densities in the diffuse (i.e. n < 1 cm" 3 ) ionized and 
atomic gas for sightlines towards ~ 200 pulsars and ~ 400 
stars within a few kiloparsecs of the Sun. The resulting 
density PDFs were consistent with a LN function over ~ 2 
dex, but the precise shape of the PDF (i.e. the dispersion and 
density corresponding to the distribution peak) varied with 
Galactic latitude. Over larger scales, neither the density nor 



column density distribution of molecular hydrogen has been 
widely investigated. In part, this is because observations of 
12 CO(7 = 1 — >• 0) emission - the most widely-used tracer 
of extragalactic molecular gas - have been limited to low 
resolution surveys of whole galaxies, or high resolution 
(often interferometric) imaging over a small fraction of a 
galactic disk. On much smaller scales, the PDF of H2 column 
density has been examined for individual molecular clouds 
in the solar neighbourhood, most ly through extinction of 
background stars (e.g. Kainulainen et al. 2009, Froebrich & 
Rowles 2010). These studies have shown that the column 
density PDFs of non-star-forming clouds are approximately 
lognormal, but that star-forming clouds exhibit power-law 
tails at high column densities, presumably due to the forma- 
tion of high density regions undergoing localised collapse. 

Our aim in this paper is to provide a quantitative description 
of the properties of the CO emission in our high angular 
resolution survey of M51's inner disk (PAWS, Schinnerer 
et al., in preparation), in order to provide simple empirical 
benchmark s for models of the H2 content of galactic disks 
(e.g. Tasker & Tan 2009 ; Dobbs et al. 201 1 ). For this purpose, 
we present PDFs of CO integrated intensity and brightness 
temperature in M51, M33 and the Large Magellanic Cloud 
(LMC), as well as CO PDFs for different dynamical environ- 
ments within the PAWS field. The extent to which the CO 
emission in M51 can be attributed to discrete, self-gravitating 
structures akin to giant molecular clouds (GMCs) in the 
Milky Way is discussed elsewhere (Colombo et al., sub- 
mitted). The rest of this paper is organized as follows. We 
summarize the origin and characteristics of the CO datasets 
that we have used in Section [2] Our method for constructing 
the PDFs is outlined in Section [3] and our results are pre- 
sented in Section [4] Our discussion in Section [5] focusses 
on the relationship between PDFs constructed from our CO 
observations on ~ 40 pc scales and the PDFs of H2 density 
and column density predicted by numerical models, and on 
the connection between GMC properties, star formation and 
the shape of the PDF for different M51 environments. We 
conclude with a summary of our key results in Section [6] We 
include three appendices, where we describe how we tested 
the robustness of the PDFs to non-physical effects such as 
resolution, sensitivity and signal identification, and assessed 
the uncertainty associated with our estimates for the slopes of 
the GMC and young stellar cluster mass distributions. 



2. DATA 

2.1. M51 

The CO data for M5 1 were obtained by the Plateau de Bure 
Arcsecond Whirlpool Survey (PAWS, Schinnerer et al., 
in preparation). PAWS observations mapped a total field 
of view of approximately 270" x 170" in the inner disk 
of M51 in the ABCD configurations of the PdBI between 
August 2009 and March 2010. Since an interferometer filters 
out low spatial frequencies, the PdBI data were combined 
with observations of CO emission in M51 obtained using 
the IRAM 30 m single-dish telescope in May 2010. In 
this section, we summarize the most important aspects of 
the observations and data reduction; the PAWS observing 
strategy, data reduction and combination procedures, and flux 
calibration are described in detail by Pety et al. (submitted). 



CO PDFs in M51 by PAWS 



2.1.1. PdBI Observations 

The PdBI observations consisted of two 30-field mosaics, 
centered such that their combination covers the inner part of 
M51. The mosaic pointings follow a hexagonal pattern, with 
each pointing being separated from its nearest neighbors by 
the primary beam full width at half maximum (FWHM). Each 
pointing was observed for 3 x 15 seconds in turn, allowing 
us to cycle completely through one mosaic pattern between 
calibrations, which were obtained every 22.5 minutes. The 
hexagonal pattern ensures Nyquist sampling along the 
declination axis but slightly undersamples the beam along the 
right ascension axis. The total telescope time in all four array 
configurations was 169 hours. The total on-source integration 
time during which useful data were obtained was 126.5 hours. 

For the PAWS observations, the two polarizations of the 
PdBI's single-sideband receiver were tuned to 115.090GHz, 
i.e. the 12 CO(7 = 1 — > 0) rest frequency redshifted to the 
LSR velocity (471.26 kms" 1 ) of M51. Four correlator bands 
of 160 MHz per polarization were concatenated to cover a 
total bandwidth of 550 MHz, corresponding to a velocity 
bandwidth of 1430 km s" 1 . The intrinsic frequency (velocity) 
channel spacing was 1.25 MHz (3.25 kms" 1 ). We later 
smoothed the data to a velocity resolution of 5 kms" 1 to 
reduce the effect of correlation between adjacent frequency 
channels and increase signal-to-noise. Inspection of the 
data showed CO emission between ±110 kms" 1 of M51's 
systemic velocity. We therefore imaged and deconvolved 
120 x 5 kms" 1 channels covering the LSR velocity range 
[174,769] kms" 1 . 

Calibration of the PdBI data was carried out using stan- 
dard methods implemented in GILDAS/CLIC^The bright 
(~ 10 Jy) quasars 0851+202 and 3C279 were used as band- 
pass calibrators. The temporal phase and amplitude gains 
were obtained from spline fits through regular measurements 
of the quasars 1418+546, 1308+326 and J1332+473. The flux 
scale was determined against the PdBI's primary flux calibra- 
tor, MWC349, and were found to be accurate within ~ 10%. 



every 12 minutes to establish the temperature scale. 

M51 observations were conducted with the EMIR receivers 
and WILMA autocorrelator backend to simultaneously record 
data for the 12 CO(7 = 1 -> 0) and 13 CO(7 = 1 -+ 0) emission 
lines. We used the upper sideband for 12 CO(7 = 1 — > 0), 
with a total bandwidth of 8 GHz. The channel spacing was 
2 MHz, corresponding to a velocity channel spacing of 5.4 
and 5.2 km s" 1 at 1 10 and 115 GHz respectively. 

The GILDAS/MIRA software was used to calibrate the 
temperature scale of the 30 m data. "OFF" spectra were 
constructed using GILDAS/MIRA's default scheme, i.e. 
averaging the closest (in time) observations together. 
These"OFF" spectra were then subtracted from the corre- 
sponding on-source spectra. Visual inspection indicated the 
presence of signal between -200 and +300 kms" 1 of M51's 
systemic velocity. A third-order polynomial was fit and 
subtracted from each spectrum. For the baseline fitting, we 
used an outlier-resistant approach and excluded regions of 
the spectrum that were known to contain bright emission, 
based on our inspection of trial data reductions or other 
observations. We experimented with higher and lower order 
baselines and found a third-degree fit to yield the best results. 
After fitting, we compared the RMS noise about the baseline 
fit in signal-free regions of each spectrum to the expected 
theoretical noise. Based on this comparison, we rejected a 
small number of spectra where the observed noise was much 
greater than expected. 

We gridded the calibrated, off-subtracted, baseline-subtracted 
spectra into a data cube with a pixel size of 4", weighting 
each spectrum by the inverse of the RMS noise. For the 
gridding, we employed a gaussian convolution kernel with a 
FWHM of 8". This gridding strategy increases the effective 
FWHM of the beamwidth to - 23 "5 at 115 GHz. After 
gridding, we fit and subtracted a second set of third-order 
polynomial baselines from each spectrum; this procedure was 
a minor refinement to the initial (pre-gridding) baseline fit. 



2.1.2. IRAM 30 m Observations 

CO emission in M51 was observed with the IRAM-30m 
single dish telescope in order to recover the low spatial 
frequency information filtered out by the PdBI. A ~ 60 
square arcminute field, covering the entire M5 1 system, was 
mapped in position-switch on-the-fly observing mode. For 
this, we divided the survey field-of-view into seven regions. 
Four regions covered the central 400" x 400" of M51; the 
remaining three regions extended the coverage to include the 
ends of the spiral arms and M51's companion, NGC5195. To 
suppress scan artifacts, each region was scanned in orthogo- 
nal directions, i.e. along the right ascension and declination 
axes, with each position in the central 400" x 400" observed 
34 times on average. The slew speed was 8"/sec and the 
dump time was 0.5s, yielding ^5.5 integrations per beam in 
the scanning direction. The scanning rows were separated by 
8", slightly oversampling the beam. We checked the pointing 
every hour, and estimate the positional accuracy to be ~ 2". 
Hot and cold loads plus the sky contribution were observed 



1 See |http: //w ww.iram. f r/IRAMFR/GILDAS for more informa- 
tion about the GILDAS software iPety 20051. 



2.1.3. Combination of Single-dish and Interferometric Data 

The final PAWS data cube is a joint deconvolution of the 
PdBI and IRAM 30 m data sets. The short-spacing visi- 
bilities not sampled by the PdBI were recreated from the 
single-dish map using the GILDAS /MAPPING software. 
For this, the map was deconvolved from the IRAM 30 m 
beam in the Fourier plane before multiplication by the 
PdBI prim ary bea m in the image plane (as described by 
Rodriguez-Fernandez et al. 2008). After a further Fourier 
transform, pseudo-visibilities were sampled between and 
15 m (the diameter of a PdBI antenna), and these visibilities 
were then merged with the interferometric observations. For 
the joint deconvolution, we used an adaption of the Htigbom 
CLEAN algorithm, as implemented in GILDAS/MAPPING. 
Supports defining the region to search for CLEAN com- 
ponents were defined for each velocity channel, based on 
where significant emission was detected in the 30 m cube. 
The convergence of the deconvolution was checked in three 
different ways. First, the cumulative flux as a function of the 
number of CLEAN components converged in each channel. 
Second, the residual channel images look like noise. Both 
criteria indicate a satisfying convergence of the deconvo- 



lution. Finally, we deconvolved the data a second time 
using exactly the same method except that we doubled the 
maximum number of CLEAN components from 320,000 to 
640,000. The subtraction of both cubes again looks like noise. 

The effective angular resolution of the final combined 
PAWS data cube is 1" 16 x 0''97, corresponding to a spatial 
resolution of ~ 40 pc at our a ssumed distance to M51 
(7.6 Mpc, Ciardullo et al. 2002). The data cube covers 
the LSR velocity range 173 to 769 kms" 1 and the width 
of each velocity channel is 5 kms" 1 . The mean RMS of 
the noise fluctuations across the survey is (Trms ~ 0.4 K 
in a 5.0 kms" 1 channel. For a typical CO linewidth of 
15 kms" 1 , this corresponds to an average sensitivity of 
3.5 Kkms" 1 for the map of CO integrated intensity. A map 
of the noise fluctuations across the PAWS field, overlaid 
with contours of CO integrated intensity, is shown in FigurefT] 



we identify significant emission using the masking method 
described by Pety et al. (submitted). For our co mpara tive 
analysis of M51, M33 and the LMC in Section |44| we 
construct an initial mask that contains pixels above a 5 a kms 
threshold over two or more contiguous velocity channels. 
This mask defines a high significance core, which is then 
expanded to include all connected pixels above \.2<jrms over 
at least two velocity channels. We discuss the rationale for 
using these blanking masks, and the influence of different 
masking techniques on the shape of the PDFs in Appendix [A] 
The total number of independent data points within the /(CO) 
PDFs varies between - 15000 and - 250000 (for the T mh 
PDFs, this increases by factor of ~ 7). For the /(CO) PDFs, 
we normalise the histogram by the number of pixels within 
the survey field-of-view, not by the number of pixels where 
significant emission is detected. Likewise, we normalise the 
r m b histogram by the number of independent (x,y, v) elements 
in the data cube. 



2.2. M33 



As we discuss in Section 4.3 the PDFs of CO integrated in- 



For M 33, we use the CO data presented by Rosolowsky et al. 
(2007), which combines observations by BIMA (Engargiola 
|et al. | |2003) and FCRAO ( [Heyer et al .12004) . The common 
field-of-view of the single-dish and interferometer surveys is 
0.25 square degrees, covering most of M33's optical disk. The 
angular resolution of the combined cube is 13"2 x 12"9, cor- 
responding to a spatial resolution of 53 pc for ou r assumed 
distance to M33 of 840 kpc (e.g. |Galletiel :al.|2004) . The data 
covers the LSR velocity range [^4-00,40] kms" 1 , and the ve- 
locity channel width is 2.0 kms" 1 . The average RMS noise 
per channel is 0.24 K. 

2.3. The Large Magellanic Cloud 

The CO data for the LMC were obtained by the Magellanic 
Mopra Assessement (MAGMA). The MAGMA survey de- 
sign, data acquisition, reduct ion procedures an d calibration 
are described in detail by [Wong et all (|20TT). MAGMA 



mapped CO cloud complexes that had been id entified at 
lower resolution by NANTEN (F ukui et al||2008) , target- 
ing 114 NANTEN GMCs with CO luminosities greater than 
7000 Kkms" 1 pc 2 , and peak integrated intensities greater 
than 1 Kkms" 1 . The combined field-of-view of the MAGMA 
survey is ^3.6 square degrees. Although the clouds tar- 
geted for mapping represent only ~ 50% of the clouds in 
the NANTEN catalog, the region surveyed by MAGMA con- 
tributes ~ 80% of the total CO flux measured by NANTEN. 
The MAGMA LMC data cube has an effective resolution of 
45", corresponding to a linear resolution of ~ 1 1 pc at the dis- 
tance of the LMC (50.1 kpc, Alves 2004). The velocity chan- 
nel width is 0.53 kms" 1 , and the total LSR velocity range of 
the cube is 200 to 305 kms" 1 . The average RMS noise per 
channel across the MAGMA field is 0.3 K. 

3. METHODS 

3.1. CO Probability Distribution Functions 

We construct PDFs of CO integrated intensity /(CO) and 
CO brightness r m t,. The /(CO) PDF is simply a histogram 
of the (x,y) pixel values within a /(CO) map, while the 
r m b PDF is a histogram of the (x,y,v) pixel values within a 
spectral line cube. The PDFs are constructed after applying 
a blanking mask that identifies genuine emission within 
the data cubes. For our analysis of the PAWS data cube, 



tensity and CO brightness T m \, within M51 exhibit diverse 
shapes that are often inconsistent with a simple functional 
form such as a lognormal or power-law distribution. We there- 
fore parameterise the shape of each PDF using the brightness 
distribution i ndex (BDI), a metric recently devised by Sawada 
et al. ( 2012b i to characterize the ratio between faint and bright 
T2 CO(7 = 1 ->■ 0) emission within an 0?8 x 0?8 field in the 
Galactic plane. More precisely, we specify the BDI of CO 
brightness as: 



BDI = log 



Er 2 



<Ti<Ti 



T< 



,Z-,r 



<Ti<T, 



T, 



(1) 



where 7] is the brightness of the i th pixel, and (7b, T\ , T 2 , T3) = 
(1.2,2.5, 5, oo)K are the thresholds that we use to define 
faint and bright emission. These are not the same thresh- 
olds adopted by |Sawada et al.| ( |2012b[ ) for their analysis at 
< 1 pc resolution, but are chosen such that variations in the 
shape of the PDF are captured (i.e. BDI is defined) for all the 
M51 environments that we analyse and that all the pixels in- 
cluded in the calculation contain significant CO emission (i.e. 
7] > 3ctrms)- We define an equivalent parameter for the /(CO) 
PDFs, which we refer to as the integrated intensity distribu- 
tion index (IDI), Analogous to Equation [Tl we specify this 
as: 

IDI = log( ^^4) (2) 



. 2/r 



</;<?! ' 



adopting (I ,I u I 2 ,h) = (10.5,25, 60,oo) Kkms" 1 . 

When appropriate, we derive the best-fitting LN function to a 
PDF using a Levenberg-Marquardt fit to the function: 



P(s) = co x exp 



-(log j -log so) 2 

2x 2 



(3) 



We define the probability P(s) as the number of pixels in the 
bin divided by the total number of pixels in the map (or cube). 
Depending on context, s represents /(CO) or T m ^. Only bins 
to the right of the peak of the PDF, brighter than 4-a RMS and 
containing ten or more pixels are used to derive the fit. Several 
of the PDFs resemble power-laws more than LN functions. In 
this case, we estimate the best-fitting slope of the power-law 
using ordinary least squares linear regression. 
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Fig. 1. — The RMS of the noise fluctuations across the PAWS survey region. The greyscale image is shown in units of T m b K. The black 
contours indicate /(CO) = 20, 70, 120, 170, 220, 320, 320 Kkm s _1 , as measured by PAWS. 



3.2. Identifying and Parameterising GMC Properties 



In Section 5.2 we investigate whether there is a connection 
between the shape of the CO PDFs and the properties of giant 
molecular clouds (GMCs) identified within the PAWS field. 
For this, we use the GMC catalog presented by Colombo 
et al. (submitted). The catalog was cons tructed using the 
CPROPS package (Rosolowsky & Leroy 2006, henceforth 
RL06). CPROPS uses a dilated mask technique to isolate 
regions of significant emission within spectral line cubes, and 
a modified watershed algorithm to assign the emission into 
individual clouds. Moments of the emission along the spatial 
and spectral axes are used to determine the size, linewidth 
and flux of the clouds. 

To generate the PAWS GMC catalog, CPROPS first identifies 
significant emission by finding pixels with CO brightness 
r m b above a Aijrms threshold across two adjacent velocity 
channels, where the RMS noise ctrms is estimated from the 
median absolute deviation (MAD) of each spectrum. This 
mask is then expanded to include all connected pixels with 
r m b > \.5o~rms- Emission regions are then decomposed into 
GMCs by identifying emission that can be uniquely associ- 
ated with local maxima. Full details of the decomposition 
procedure are presented in Colombo et al. (submitted). 

We adopt the default CPROPS definitions of GMC proper- 
ties. The cloud radius is defined as R = 1.91oRpc, where ctr 
is the geometric mean of the second moments of the emis- 
sion along the cloud's major and minor axes. The velocity 
dispersion er v is the second moment of the emission distribu- 
tion along the velocity axis, which for a Gaussian line profile 
is related to the FWHM linewidth, Av, by Av = V81n2cr v . 
The CO luminosity of the cloud Leo is the emission inside 



the cloud integrated over position and velocity, i.e. 

Leo [ Kkm s" 1 pc 2 ] = D 2 ( ) ZT6v6x8y , (4) 

F V 180x3600/ 

where D is the distance to the galaxy in parsecs, 8x and 8y are 
the spatial dimensions of a pixel in arcseconds, and Sv is the 
width of one channel in km s" 1 . The mass of molecular gas 
estimated from the GMCs CO luminosity Mqo is calculated 



Mco [M ] = 4.4 ^° — r L C0 

2x lO^tcm-^Kkms" 1 ) ] 



(5) 



where Xqo is the assumed CO-to-H2 conversion factor, and a 
factor of 1.36 is applied to account for the mass contribution 
of helium. The fiducial value of Xco used by CPROPS 
is Xqo = 2.0 x 10 2() cm" 2 (Kkm s" 1 )" 1 . The virial mass is 
estimated as M v ; r [M ] = 1040ct 2 /?, which assumes that 
molecular clouds are spher ical with truncated p ex r~ l density 
profiles ( |MacLaren et al.|[l988] >. CPROPS estimates the 
error associated with a cloud property measurement using a 
bootstrapping method, which is described in section 2.5 of 
RL06. 

The final PAWS GMC catalog contains 1507 objects. The 
GMCs have peak brightness temperatures between ~ 2 and 
16K, radii between 5 and 150 pc, and velocity dispersions 
between 1 and 30 kms" 1 . The catalog and the properties of 
GMCs in different environments within the PAWS field are 
the subject of a companion paper (Colombo et al., submitted). 
In general, the physical properties of the cataloged GMCs 
are similar to the GMCs identified by CO surveys of the in- 
ner Milky Way and other nearby galaxies, although GMCs in 
M51 tend to be larger, brighter and have higher velocity dis- 
persions and mass surface densities relative to their size than 



the GMCs in nearby low-mass systems such as the Magellanic 
Clouds and M33 (Hughes et al., in preparation). The spa- 
tial resolution and sensitivity of PAWS is sufficient to resolve 
structures with size a nd mass co mparable to a typical Galactic 
GMC (50pc, 10 5 M |Blitz|l993} . We note, however, that CO 



emission is almost ubiquitous across the PAWS field, and that 
much of the emission resides in large (^kpc-sized) regions of 
high brightness that bear little resemblance to Galactic GMCs. 
Overall, the cataloged GMCs account for approximately half 
of the total CO flux within the PAWS datacube, a fraction that 
varies from ~ 40% in the interarm region to ~ 60% in the 
spiral arms and central zone. 

4. RESULTS 

4.1. I(CO)PDF 

The PDF of CO integrated intensity for the entire PAWS 
field is presented in Figure [2ta]. The distribution is ad- 
equately described by a LN function, with a mean of 
(/(CO)) = 17.4 Kkms" 1 and a logarithmic width of 0.44. 
The logarithmic dispersion in the fit residuals e for bins above 
the 3cr sensitivity limit that contain more than five counts is 
0.08. Relative to this LN function, there is some evidence 
for a truncation at high /(CO) values (> 200 Kkms" 1 ). In 
principle, this could be due to opacity of the 12 CO(7 =1—5-0) 
emission line, an effect that is often observed on parsec 
scales in regions of high extinction (Ay > 5- 10 mag, e.g. 



Lombardi et al. 2006; Pineda et al. 2008 ), but has rarely been 
considered for the scales probed by extragalactic observations 
(cf. Dickman et al.|l98 6). Against this intepretation, the 99 th 
percentile of the CO peak brightness within the PAWS field 
is only ~ 7 K, suggesting that there are negligible sightlines 
where the CO emission completely fills the telescope beam. 
Furthermore, recent numerical simulations that examine the 
ability of /(CO) to trace the H2 column density on 20 to 60 pc 
scales show that saturation tends to produce a secondary peak 
in the /(CO ) PDF rather than a sm ooth truncation ( S hetty 
|et al.||2(7TTj |Feldmann et~aT] |2"0"T2] see also Section [54}. 
The finite resolution of observational data can also produce 
a truncation at high intensities (we explore this effect in 
Appendix [Bj, In this case, however, we would expect the 
/(CO) PDFs for subregions within the PAWS field to exhibit 
similar thresholds, whereas severa l of them are consistent 
with pure LN functions (see SectionpOjl. 



The truncati on of the PDF is Figure [2fa] may therefore 
be physical. Elmeg reeii] ( |2011]l show that the density PDF 
should fall beneath a pure LN function at high gas densities 
if the Mach number decreases with increasing average gas 
density (as would be expected for a cloud that obeys the 
Larson (|1981|l size-linewidth relation). Alternatively, the 



truncation may reflect the efficacy of feedback processes 
that prevent the molecular gas from reaching very high 
mass surface densities (Sh, > 4OOM pc~ 2 , assuming 
X C o = 2.0 x 10 20 cm" 2 (Kkms -1 )" 1 ). We discuss physical 
processes that could be influencing t he sh ape of the /(CO) 
PDFs in M51 in more detail in Section [ 



4.2. T mb PDF 

The PDF of CO brightness for the PAWS cube is shown 
in Figure [2Tb]. The r m b PDF is less like a LN function 
than the /(CO) PDF, with two roughly flat segments across 



1 < T mb < 3 K and T mb > 5 K. Our LN fit to the PDF yields 
a mean (r m b) ~ 0.9 K and logarithmic width x ~ 0.35. 
Assuming that the true distribution of T mb values is LN, there 
are fewer high brightness pixels than would be expected from 
this LN function. The truncation begins to occur at a CO 
brightness temperature of ~ 5 K, which would seem too low 
to be due to opacity effects. Instead of an LN function, a 
broken power-law with a slope of ~ -1.4 for 1 < T mb < 5 K 
and a much steeper slope of ~ -4.9 for T mb > 5 K - or, alter- 
natively, a pure power-law with a truncation at ~ 5 K - may 
provide a better description of the PDF. The fit parameters 
and goodness-of-fit for the best-fitting LN and power-law 
functions to the T mb PDF in Figure [2jb] are listed in Table [T] 



4.3. M51 Environments 

An important question that we would like to address with 
the PAWS data is whether the organization and physical 
properties of molecular gas depend on galactic environment. 
Variations in the shape of the PDF could reflect differences 
in the relative importance of self-gravity, star formation 
feedback or gas flow in different parts of the galactic disk, 
which in turn might influence the ability of the molecular gas 
to form stars. Within the PAWS field, there are three main 
regions where the gas is likely to experience distinct physical 
conditions: within the strong spiral arms, the inter-arm region 
situated upstream and downstream of the spiral arms, and the 
central region, where the gas is influenced by the presence of 
a nuclear stellar bar (Zari tsky et al.|1993] >. These regions can 
be further classified according to their level of star formation 
activity (as traced by e.g. H a) and/or gas flows, which we 
determine using the present-day torque profile (Meidt et al., 
submitted). 

Here, we analyse seven regions within the PAWS field where 
we expect the molecular gas to experience different dynam- 
ical effects (see Figure B). We define the different spiral 
arm regions according to the direction of gas flows driven 
in response to the underlying gravitational potential, which 
we derive from a map of M51's stellar mass distribution 
(Meidt ~et al.|2012 i. The width of the spiral arms are defined 
with respect to observed gas kinematics. We determine the 
zone of enhanced spiral streaming centered around the arm 
by measuring the (rotational) auto-correlation of azimuthal 
streaming velocities in the PAWS field (Colombo et al., 
submitted). We construct azimuthal profiles of the auto corre- 
lation signal in a series of radial bins and take the width of the 
signal at 95% maximum as our measure of the kinematic arm 
width. The average kinematic width from along the two arms 
is centered on the spiral arm ridge line, defined by eye using 
the PAWS map of CO peak brightness. Both the location of 
the ridge and the width are assumed to be symmetric. The 
interarm region is divided into upstream and downstream by 
the midpoint of the spiral arm ridge lines. The definition of 
the spiral arm regions is based on the identificatio n of distinct 
spiral patte rns withi n the galactic disk (cf. Vogel et al. 1993, 
Shetty et al.p07]|Meidt et al |2008|pobbs et al |2010| > that 



we refined and describe in detail elsewhere (Meidt et al., 
submitted; Colombo et al., submitted). 

The seven zones that we use to conduct our analysis are: 

• Nuclear bar: the region at galactocentric radii R < 23". 
The boundary is defined by the bar corotation reso- 



CO PDFs in M51 by PAWS 



S = log(E H 2/[M S ol PC" 2 ]) 

1.0 1.5 2.0 2.5 3.0 



s = log(E H 2/[M so i per 2 ]) 

0.5 1.0 1.5 2.0 




0.5 1.0 1.5 2.0 2.5 -0.5 0.0 0.5 1.0 

s = log(Ico/[K km s" 1 ]) s = log(T mb /[K]) 

Fig. 2. — The [a] /(CO) and [b] T m b PDFs for emission within the PAWS field. In both panels, a dashed parabola indicates the LN function 
that provides the best-fit to the PDF. In panel [b], the dot-dot-dashed lines represent two segments of a broken power-law, which also provides a 
reasonable fit to the distribution. The grey shaded region indicates values beneath our nominal 3o-rms sensitivity limits of 10.5 Kkm s - ' (panel 
[a]) and 1.2 K (panel [b]). The top horizontal axis shows the equivalent H2 mass surface density for the /(CO) or T m b value on the lower 



axis, assuming Xco : 



associated with simple counting (yN) errors. 

nance, inside of which the bar exerts negative torques 
and drives gas radially inwards. 

• Molecular ring: the region 23 < R < 35". Here the 
gas is influenced by both the bar and innermost spi- 
ral arms. Outside the bar corotation resonance, gas is 
driven radially outwards, while the spiral drives gas ra- 
dially inwards inside its own corotation. These oppos- 
ing torques accumulate gas in a ring-like structure. The 
region hosts some of the most active high-mass star for- 
mation in M5 1 . 

• Inner density-wave spiral arm: the arm region 35 < 
R < 55". The inner boundary is defined by the molec- 
ular ring, and the outer boundary is the corotation ra- 
dius of the density-wave spiral arms. Within this zone, 
gas is driven radially inward by negative spiral arm 
torquing. Despite the high gas surface densities in this 
region, there is little star formation - as traced by H a 
and 24 /im emission - that is directly associated with 
the brightest CO emission (Schinnerer et al., in prepa- 
ration). 

• Outer density-wave spiral arm: the arm region 55 < 
R < 85". This region extends from the density-wave 
corotation resonance to the start of the material spiral. 
Within this zone, gas is driven radially outward by pos- 
itive spiral arm torquing. 

• Material spiral arm: the arm region R > 85". This re- 
gion extends from the boundary of positive arm torques 
associated with the density wave spiral to the edge of 
the PAWS field. There is some indication that gas flows 
radially inwards in this zone. 

• Interarm region, downstream of the spiral arms. 

• Interarm region, upstream of the spiral arms. 



' (Kkms ') ' and a helium contribution of 1.36 by mass. The vertical error bars represent the uncertainty 



Finally, we note that the projected area of the seven regions is 
still quite large (between ~ 2 and 17kpc 2 ) and each contains 
a statistically significant number of GMCs (> 100). The 
PDFs of CO emission in these regions are therefore more 
comparable to the PDFs of simulated galactic disks than the 
PDFs of individual clouds. 

Figure [4] shows that the /(CO) PDFs for different M51 
environments exhibit diverse shapes. The panels of Figure |4] 
are ordered such that the PDF amplitudes decrease from top 
left to bottom right, which reflects the fact that CO emission 
is more prevalent in the arms and central region of M51 
than in the interarm region. The PDFs also tend to decrease 
in width, indicating that the fraction of pixels with bright 
CO emission declines with the overall frequency of CO 
detections. The PDFs of the spiral arms are notably wider 
than for the interarm environments; it is also evident that the 
PDF corresponding to the first spiral pattern (Al, i.e. the 
density wave spiral arm) is wider than the PDF for the second 
spiral (A2, the material arm). Since /(CO) is the integral of 
the CO brightness over the line profile, this variation in the 
PDF width would seem consistent with the results of previous 
studies (e.glGarc ia-Buril lo et al.|19 93 1 Kuno & Nak ai|1997| 
|Aalto et al.||1999[ |Schuster et al.||2007| > that find that the 
average CO integrated intensity and CO linewidth decreases 
with increasing distance along the arms, and from the arm to 
the interarm region. 

The differences in the width of the PDF between M51 
environments are reflected in the IDI values that we derive, 
which become more positive as the number of pixels with 
/(CO) > 60 Kkms -1 increases (see Table J2|i. The devel- 
opment of more high brightness emission appears to be 
accompanied by a change in the PDF shape: a LN function is 
a better description of the PDFs in the interarm region than 
in the spiral arm and ring regions. The /(CO) PDFs for the 
first spiral pattern, especially inside corotation (All), appear 




202.58 202.5 202.48 202.46 202.44 202.42 
RA (degrees) 



Fig. 3. — The different environments within the PAWS field that we analyse in this paper, indicated using grey shading in each panel, (a) 
Nuclear bar (bar); (b) Molecular ring (ring); (c) First spiral pattern inside corotation (All); (d) First spiral pattern outside corotation (AlO); 
(e) First spiral pattern (density wave arm, Al); (f) Second spiral pattern (material arm, A2); (g) Interarm region upstream of the spiral arms 
(up); (h) Interarm region downstream of the spiral arms (down). The black contours in all panels indicate /(CO) = 25 Kkms~', as measured by 
PAWS. 
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more like broken or truncated power-laws than LN functions. 
The PDFs in the center of M5 1 also diverge from a LN shape: 
the /(CO) distribution in the molecular ring is essentially flat 
between ~ 20 and 150 Kkm s" 1 , while the PDF of the bar is 
the only region with an unambiguous decline at low intensi- 
ties (/(CO) < 50 Kkms" 1 ). Several PDFs appear truncated 
near /(CO) ~ 300 Kkms" 1 : this is seen most clearly in the 
molecular ring, but the distributions in the first spiral arm 
regions also decline steeply for /(CO) > 300 Kkms -1 . 

In Figure [5] we present the PDFs of CO brightness for the 
different M51 environments. The distributions are more 
uniform than those of integrated intensity, but variations 
similar to those identified for the /(CO) PDFs are still 
evident. The PDF amplitude tends to decrease from panel 
[a] to [h], and only the interarm regions and second spiral 
arm (A2) yield PDFs that are approximately LN across the 
observed range of r m b values (see Table [TJ. As noted for the 
/(CO) PDFs, regions with a relatively wide T mb PDF (e.g. the 
ring, bar and spiral arms) have more positive BDIs and also 
tend to diverge from a LN shape, in this case developing a 
pronounced change of slope near T m b > 6 K. This effect is 
most clearly seen for the PDFs of the molecular ring (panel 
[b]) and the first spiral pattern inside corotation (panel [c]), 
but generally it appears that an increase in the fraction of 
high brightness CO emission is associated with a PDF that 
more resembles a truncated power law than a LN function. 
We discuss this result in relation to similar trends observed 



for PDFs of CO emission in the Galaxy in Section 5.1 



4.4. Comparison between M51, M33 and the LMC 

Finally, we can compare the PDFs of CO integrated intensity 
and CO brightness for the inner disk of M51 to the corre- 
sponding PDFs for other nearby galaxies. As we discuss 
in Appendix IB] the shape of the PDF is sensitive to the 
resolution andsensitivity of the data. Prior to constructing the 
PDFs, we therefore degraded the M51 and LMC data cubes 
to the same spatial resolution as the M33 cube (~ 53 pc), and 
folded the M33 and LMC data cubes along the velocity axis 
to the same channel width as the M51 cube (5 kms~ ). We 
interpolated all the cubes onto an (x,y) grid with the same 
pixel dimensions in physical space (15 x 15 pc). Significant 
emission was identified according to the method outlined 
in Section [3] The resulting masks were applied to the 
original data cubes, and the integrated intensity images were 
constructed by summing unblanked pixels across the full 
velocity bandwidth of each survey. 

The PDFs obtained from the /(CO) maps of M33, the LMC 
and M51 are shown in the left panels of Figure [6] The shape 
of the /(CO) PDF for M33 is highly uncertain due to the 
modest sensitivity of the BIMA+FCRAO data cube. Never 
the less, it is obvious that CO emission in the LMC and M33 
are alike in the sense that the maximum observed /(CO) 
intensities are ~ 10 Kkms -1 at the resolution of our analysis, 
and not a few times 100 Kkms" 1 , as observed for M51. The 
/(CO) PDF for the LMC appears to be well-represented by 
a narrow LN function with mean (/(CO)) = 2 Kkms" 1 , and 
a logarithmic width x = 0.3. Since MAGMA is a targeted 
rather than spatially complete survey of the LMC disk (see 
Section EJ, the amplitude of its PDF is biased high compared 
to that of M51 and M33; normalizing by the full projected 



area of the LMC's Hi disk, rather than the MAGMA field- 
of-view would reduce the amplitude by more than an order 
of magnitude. The MAGMA survey strategy of targeting 
the brightest clouds in the NANTEN catalog also means that 
the shape of the PDF is biased towards high CO intensities; 
extending MAGMA to fainter clouds would recover a greater 
fraction of pixels with low CO brightness and narrow the 
PDFs in Figure [6jc] and [f]. Relative to observed /(CO) 
PDF of the LMC, the M51 /(CO) distribution peaks at higher 
CO intensity, (/(CO)) ~ 12 Kkms" 1 , and is also wider by a 
factor of ~ 2 in the logarithm (see Table [T}. The difference 
between the best-fitting LN function derived for the M51 
distribution in Figure[6ta] and that in Figure 2ta] is consistent 
with the differences that we observe for different masking 
techniques for identifying significant emission within the 
data cube (see Appendix |A). 

The PDFs of CO brightness for the three galaxies are pre- 
sented in the right panels of Figure [6] For M5 1 , we fit the 
distribution of CO brightness with a LN function with mean 
(Zmb) = 0.8 K, logarithmic width x = 0.4. This is consistent 
with the best-fitting LN function derived for the PAWS data 
in Figure EJb], which used a more sophisticated masking 
technique to identify significant emission. Alternatively, a 
broken power-law with a shallow slope of ~ —1.3 between 
r m b = 1 and 4 K, and a much steeper slope (~ — 7. 1) above 4 K 
also fits the T m \, distribution for M51 reasonably well. For the 
LMC, the best-fitting LN function has mean (T mb ) = 0.1 K 
and logarithmic width x = 0.3; there is no sign of a truncation. 
A simple power law with a slope of ~ -2.2 also adequately 
represents the distribution. There are insufficient pixels 
with significant emission in the M33 data cube to attempt 
to fit the r m b PDF with a LN function. A simple power-law 
with a slope of ~ -3.6 provides a reasonable fit to the 
distribution for pixel values r m b > 0.3 K. We note that the 
slopes of the power-laws that fit the T m \, PDFs reproduce the 
trends observed for the giant molecular cloud (GMC) mass 
distribution in the three galaxies, i.e. a shallow slope for M51 
(below r m b ~ 4 K) and a much steeper brightness distribution 
for the low-mass galaxies (Hughes et al., in preparation). We 
discuss the connection between the shape of th e CO PDFs 
and the GMC mass distribution further in Section [ 



5. DISCUSSION 

5.1. Comparison to Previous Observations 

Both high resolution and wide-field coverage are necessary to 
characterise the CO emission in galaxies on spatial scales that 
are relevant for star formation, hence few extragalactic stud- 
ies have produced PDFs of CO brightness and/or integrated 
intensity that represent a significant fraction of a galactic 
disk. O ne exception is an analysis of the LMC by Wong et al. 
( |2"0TT] >, which found that on - lOpc scales the /(CO) PDF 
was roughly consistent with a narrow lognormal function 
((E#,) = 16M Q pc" 2 , x ~ 0.3 dex) at high column densities. 
The authors noted some evidence for a truncation around 
Eh, = 200 M Q pc" 2 , which they tentatively attributed to opac- 
ity effects in the 12 CO(7 = 1 — > 0) line. Notably, however, the 
MAGMA data show no evidence for a power-law excess at 
high column densities, as has been observed on parsec scales 



within star-forming Galactic clouds (e.g. Kainulainen et al. 



2009 ), suggesting that the structure of LMC molecular clouds 
is still dominated by turbulence on ~ 10 pc scales. 
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Fig. 4. — The /(CO) PDFs for different regions within the PAWS field. The grey shaded region represents values beneath our nominal 3a RM s 
sensitivity limit of 10.5 Kkms~'. Where a LN (power-law) function provides an good description of the PDF, it is indicated by a dashed 
(dot-dot-dashed) line. The vertical error bars represent the uncertainty associated with simple counting (v^V) errors. 
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Fig. 5. — The PDFs of CO brightness for different regions within the PAWS field. In each panel, the grey shaded region represents values 
beneath our nominal 3a sensitivity limit of 1.2 K. When appropriate, the best-fitting LN (power-law) function is indicated by a dashed (dot-dot- 
dashed) line. The vertical error bars represent the uncertainty associated with simple counting (\/A0 errors. 
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TABLE 2 

Brightness and Integrated Intensity Distribution Index for M51 Environments 



Region 
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BDI 


IDI 




[10 7 Kkms _1 pc 


2 ] 




Global 


70.4 


-0.66 


0.40 


Nuclear Bar 


6.6 


-0.85 


0.76 


Molecular Ring 


16.6 


-0.13 


1.08 


Arm 1 inside corotation (All) 


11.9 


-0.59 


0.52 


Arm 1 outside corotation (AlO) 


17.6 


-0.83 


0.55 


Arm 1 (Al) 


29.6 


-0.72 


0.54 


Arm 2 (A2) 


6.7 


-0.76 


0.12 


Upstream 


4.7 


-1.73 


-0.95 


Downstream 


6.2 


-1.50 


-0.65 



The total CO luminosity (column 2), brightness distribution index 
(BDI, column 3) and integrated intensity distribution index (IDI, 
column 4) for the different M51 environments (see Figurepl. The 
BDI and IDI values are calculated according to Equations^and[2] 
respectively. More positive BDI and IDI values indicate PDFs that 
have a larger fraction of pixels at high CO intensities. 
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Fig. 6. — PDFs of /(CO) (panels [a] to [c]) and T m t (panels [d] to [f]) for M51, the LMC and M33, constructed using a dilated mask technique. 
The data sets have been smoothed to the same spatial scale, and interpolated onto an (x,y, v) grid with the same spatial and spectral dimensions. 
In all panels, the grey shaded region corresponds to pixels beneath the Ixjrms sensitivity limits of the individual data cubes. The vertical 
error bars represent the uncertainty associated with simple counting {\/N) errors. Where relevant, the best-fitting LN (power-law) function is 
indicated by a dashed (dot-dot-dashed) line. 



The CO emission in M51 itself has be e n analysed by nu- 
merous authors (e.g. |Vogel et al.]|1988| |Rand||1993| |Aalto 
etaLlT999lpelfer et al. 2003; |Shetty et aTl2 007). WitnTEe 
exception of the recent survey by |Koda et al.| (|2009), how- 
ever, observations with high spatial resolution have mostly 
focussed on a spiral a rm segment (e.g. Schinnerer et al. 
2010 Egusa et al.|201 1 1, while studies covering a significant 



et al.| [2007) have had a resolution of a few hundred parsecs 
or greater, i.e. insufficient resolution to resolve individual 
GMCs. Rather than PDFs, these lower-resolution studies have 
typically examined radial trends in properties such as the gas 
velocity dispersion, Toomre's Q parameter and the molecular 
gas depletion time r H , (e.g. Schuster et al. 2007 ; Hitschfeld 
et al.| |2009'l. From "the PDFs in Section |4.3f we would 
expect to observe a radial decline in the average value of 
fraction of the disk (e.g. Garcia-Burillo et al.|1993 Schuster Eh^ constructed from azimuthal averages, since the fraction 



14 



of bright CO emission (as parameterised by the IDI values) 
decreases along the spiral arms and also because the interarm 
region occupies an increasing fraction of the disk area with 
increasing galactocentric radius. Since the mass surface den- 
sity is an important input for the determination of Toomre's Q 
and th 2 , radial trends in these quantities may likewise reflect 
a combination of differences between the arm and interarm 
zones and variations along the spiral arms (an interpretation 
that would seem to be supported by the map o f Toomre's Q 
presented in figure 15 of Hitsc hfeld et al.| | |2009) for example). 

More generally, we note that differences in the basic prop- 
erties of the CO emission (i.e. peak brightness, velocity 
dispersion) between the arm and interarm regions of M51 
have been reported by several previous studies, and often 
interpreted as evidence for changes in the physical state 
of the molecular gas as it passes through the spiral arms. 
While the precise identification of M5 1 's arm and interarm 
zones varies (usually because M51's gaseous spiral arms 
appear wider at lower spatial resolution), CO emission in the 
interarm has been shown to have lower velocity dispersion 



etal. 



(e.g.|G arcia-Bur illo et al |1993||Aalto eTaTT"l 999; Hit schfeld 



2009), lower peak brightness (e.g. Garcia-Burillo et al. 

Tosaki et al. 2002 1, lower star formation efficiency 

(as inferred from the ratio of Ha to n CO(J = l—>0) 



1993 



emission, e.g. Rand] [T993] Tosaki et al. 20021, higher 
12 CO(7 = 1 -)■ 0)/ 13 CO(7 = 1 -> 0) isotopic ratios (e.g. |Tosaki 
et al.||2002| > and lowe r 12 CO(7 = 2 -> l)/ 12 CO(7 = 1 ->• 0) 
transitional ratios (e.g. Kodaet al.|2012 i than emission in the 
spiral arms. Most of these results suggest that molecular gas 
in the interarm region has a lower characteristic density than 
gas within the spiral arms. 

Spatial resolution is rarely a limitation for studies of molec- 
ular gas in the Milky Way, although previous analyses of 
Galactic CO emission ha ve tended to focus on the physical 
properties o f GMCs (e.g. Solomon et al. 1987| |Roman-DuvaI| 
|et al.| |2010| i, which were quickly recognized to be the 
preferred - perhaps only - site of high-mass star formation 
in the Galaxy. Recent work has emphasized, however, that 
faint spatially extended CO emission contributes s ignificantly 
to a region's total CO flux (e.g. |Goldsmith et al.||2008 



Heyer et al . 2009, L iszt et al.|2010) . Very recently, |Sawa3a1 
et al.| (|2012bt have presented PDFs of 12 CO(7 = 1 -> 0) 



and "CO(/ = 1 -> 0) brightness for an 0?8 x 0?8 field 
towards the Galactic plane at / ss 38°, observed using the 
Nobeyama Radio Observatory (NRO) 45 m telescope. These 
authors find clear differences between the PDFs constructed 
from the emission at radial velocities corresponding to the 
Sagittarius arm and those corresponding to the interarm 
regions, showing that the structural proprties of the molecular 
gas vary in response to Galactic structure. They conclude 
that compact, high brightness CO structures develop down- 
stream of the molecular spiral arms, where they are spatially 
coincide nt with signatures of active star formation (e.g. Hn 
regions). Sawada et al. (2012b) point out that their result is a 
rediscovery of a conclusion that had already been drawn by 
earlier studies. Sanders et al. (1985), for example, observed a 
connection between the location of Galactic molecular clouds 
in longitude-velocity space and their peak brightness temper- 
ature: "hot" (i.e. high brightness) clouds were preferentially 



similar scenario for a ~ 2 kpc segment of M5 1 's inner spiral 
arm, showing that both high mass (~ 1O 6 M ) CO clumps 
and H n regions are preferentially located downstream of the 
spiral arm ridge line (note, however, that the high brightness 
CO structures described by Sawa da et al.| ( |2012b] l occur on 
much smalle r spatial scales than the structures observed by 
Egusa et al.|p0TT]> in M51). 



Although the spatial resolution of the PAWS data is con- 
siderably worse than that of Galactic surveys, the PDFs in 
Figure |5] show similar trends as those repo rted by |Sawada| 
et al. (2012b). As we noted in Section |4~3*] the PDFs of the 
interarm region resemble narrow lognormal functions, while 
the PDFs in the central and spiral arms regions reach higher 
maximum intensities and tend to be better represented by 
broken power-laws. These variations in shape are reflected by 
the BDI and IDI values: low brightness emission dominates 
the total flux in both arm and interarm environments, but the 
relative contribution from bright emission increases in the 
spiral arms. Bright emission is most dominant in the center 
of M51, where ^25% of the total CO flux arises from pixels 
with r m b > 4K (by comparison, less than 5% of the emission 
in the interar m is brighter than 4K). Similar to |Sawada| 
|et al.| ( |2012b| l, we find that the BDI and IDI values are 
higher on the downstream side of the spiral arms than on the 
upstream side. Tracers of high-mass star formation, e.g. H a, 
24 /im and far ultra-violet (FUV) emission, also appear to be 
preferentially located downstream of arms (Schinnerer et al., 
in preparation), again consistent with the Galactic results. We 
discuss the connection between star formation and the shape 
of the CO PDFs for different M51 environments in more 
detail in SectionT: 



located in the sp ir al arms traced by H n regions ( George lin & 
Geor gelin||1976| l. Egusa et al. ( 2011| i present a qualitatively 



The fact that a similar relationship between CO emission 
properties and spiral arm structure is observed in both M51 
and the Milky Way w ould seem to support the ar gument by 
|Sawada et all ( |2012b) and |Sawada et aT] P012a| ) that the 
arm-interarm variations they observe reflect genuine changes 
in the density distribution of the Galactic CO-emitting gas. 
One important caveat, however, is that the PAWS data has 
much lower spatial resolution (~ 40 pc) than the Galactic 
NRO data (< 1 pc). In particular, the lower CO brightness 
temperatures that we observe in M51 (r m b < 16 K) indicate 
that our PAWS measurements reflect a combination of the 
average kinetic temperature and the filling factor of the 
CO-emitting gas within a resolution element. Variations in 
CO brightness on < 1 pc scales, by contrast, should mostly 
track variations in gas temperature and/or density since beam 
dilution should be minimal on these scales. Some of the 
variation between high and low BDI values in M51 will 
reflect changes in the covering fraction of the CO emission 
for the different M51 environments, as well as differences 
in the intrinsic brightness temperature of the CO-emitting 
structures that are m ore directly co mparable to the variations 
described by |Sawada et al.| ( [2012b| l. 

5.2. Comparison with GMC properties 

In Sections |4.3| we described variations in the characteristic 
shape of the CO PDFs for different M51 environments. To 
what extent are these differences manifested in variations 
of the properties of GMCs within each environment, or of 
the ensemble properties of a GMC population (e.g. its mass 
distribution)? Intuitively we would expect some connection 
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between GMCs and the presence of high brightness CO 
emission, since most methods for identifying GMCs from 
CO data cubes invoke either a brightness threshold or local 
maximum in the CO brightness distribution in order to define 
cloud structure. The connection may be rather indirect, 
however, since the fraction of CO emission above the PAWS 
sensitivity limit that is associated with the observationally- 
defined GMCs varies between 40 and 65%, depending on 
galactic environment (Colombo et al., submitted). 

We examined the relationship between GMC properties and 
the shape of the CO PDFs using the cloud catalog presented 
by Colombo et al. ( subm itted) and the BDI and IDI values 
calculated in Section [43] Since the BDI and IDI themselves 
exhibit a tight one-to-one correlation, for simplicity we 
refer only to the IDI in the following sections. We illustrate 
some of these correlations in the left column of Figure [7] 
Environments where bright CO emission is more dominant 
(i.e. with more positive IDI values) are associated with a 
higher maximum GMC mass M gmc $s, which we estimate 
using the 95th percentile of the GMC virial mass distribution 
(panel [a]), a greater number surface density of GMCs 
Afg mc (panel [b]), and a higher average surface density for 
individual GMCs (E Ho ) (panel [c]). The IDI is also strongly 
correlated with the slope of the GMC mass spectrum j gmc : in 
environments with more bright emission, the mass spectrum 
is shallower (panel [d]). For observations with low resolution 
(i.e. where a single resolution element is much larger than 
the characteristic size of a GMC), a good correlation between 
the prevalence of bright CO emission and the mass and mass 
surface density of identified cloud structures might arise sim- 
ply due to higher filling factors of CO emission, i.e. increases 
in the measured CO integrated intensity reflect a greater 
number of CO-emitting clouds within the telescope beam, 
rather than changes in the intrinsic properties of GMCs. We 
do not consider this to be the cause of the good correlations 
in Figure 17] however, since the PAWS resolution (~ 40 pc) 
is well-matched to the chara cteristic siz e of an individual 
Galactic GMC (50 pc, e.g. | Bhtz|[l993| ), and considerably 
less than the typical spacing between the identified GMCs (a 
few times 100 pc or greater, Colombo et al. submitted). The 
peak CO brightness temperatures of the GMCs range from 
~ 2 to 16 K, which is comparable to the values o bserved for 
Galactic GMCs (5 to 10K, |Solomon et al.|[T987l >. Since the 
molecular gas in M5 1 GMCs appears to hav e a similar kinetic 
temperature as in Galactic GMCs (~ 10 K, Schinnerer et al.| 
2010), this again suggests that the filling factor of the CO 
emission in M5 1 GMCs is close to unity. 

It is remarkable that the GMC properties are often more 
strongly correlated with the shape of the CO PDFs than 
other quantities with which they might also be expected to 
correlate. In particular, we note that M gmc .95 and j gmc are 
more tightly correlated with the IDI than with the total CO 
luminosity or the total number of GMCs in each region (the 
latter plots are not shown). This would seem to confirm 
that an increase in the maximum GMC mass is not simply 
due to an increase in the available gas reservoir and more 
adequate sampling of the top-end of the GMC mass function 
(i.e. a size-of-sample effect), and that the good correlation 
between the IDI, M gmc ^s and 7 gme arises because the density 
distribution of the molecular ISM plays a role in regulating 
the GMC mass distribution. It is also noteworthy that the IDI 
increases with both the number density of GMCs N gmc and 



the average surface density of the individual clouds (Ehi)- 
This suggests that a distinction that is sometimes drawn by 
empirical studies of extragalactic star formation between an 
increase in the number of GMCs per resolution element and 
variations in the Fb surface density on the scale of individual 
clouds (e.g. Bigiel et al. 2008) is somewhat artificial: at least 
in the inner disk of M51, clouds in environments with more 
GMCs per unit area also tend to have higher average surface 
densities. 

The good correlation between the IDI and 7 ?mc , in panel [d] 
of Figure [7J is especially noteworthy. Considerable theoreti- 
cal and observational effort has been devoted to showing how 
the shape of the stellar initial mass function might be inher- 
ited from the density structure of interstellar gas (e.g. Hopkins 
2012, Chabrier & Hennebelle 2010 and references therein), 
with many studies adopting the shape of the GMC mass func- 
tion as a description of the latter. A major problem with us- 
ing GMC mass spectra for this purpose, however, is that the 
decomposition algorithm has a major impact on the identi- 
fication and parameterization of cloud structures and hence 
the shape of the resulting mass distribution (e.g. Won g et ah] 
|2011|[Reidetal. 2010). Moreover, many widely-used decom- 
position methods are not flux conservative, discarding a con- 
siderable fraction of the CO emission that is unambiguously 
detected within a spectral line data cube. As a description of 
how dense gas is distributed within galaxies, PDFs avoid these 
ambiguities even though, as we show in the Appendix [5] the 
resolution of the observational data must be well-matched to 
the physical scales of interest in order to capture the shape of 
the PDF accurately. The PDF also conveys no information 
about the characteristic size of dense gas structures, more- 
over, so a more complete description of the organization of 
the dense ISM strictly requires an analysis of the CO PDF in 
conjunction with a metric such as the spatial power spectrum 
(as suggested by e.g. [Bournaud et al. 2010). Nonetheless, 
the plots in the left column of Figure 17] indicate that there is 
a strong relationship between the shape of the PDF and the 
mass distribution and properties of GMCs within M51 envi- 
ronments, suggesting that in real galactic disks the presence 
of bright emission and the development of massive molecu- 
lar structures are physically linked. Testing whether a sim- 
ilar connection between the GMC mass function and shape 
of the PDF holds across a range of galaxy types is a project 
that should become feasible once ALMA acquires cloud-scale 
imaging of the CO emission across the full galactic disk for a 
large sample of nearby galaxies. 

5.3. Comparison with properties of stellar clusters 

Another motivation for constructing the CO PDFs across 
a range of M51 environments is to assess whether there 
are connections between empirical tracers of star formation 
activity and the density distribution of molecular gas. Several 
empirical calibrations for the star formation (SF) rate exist 
in the literature (for a detailed comparison of the limitations 
and assumptions of different methods, see Leroy et aLJ2 012), 
but here we restrict our analysis to comparing the shape 
of the CO PD Fs to the p roperties of young stellar clusters 
identified by Chandar et al. (201 1 ) using multi-colour images 
of M51 obtained by the Advanced Ca mera for Surveys on 
board the Hubble Space Telescope (Mutc hler et al.| [2005). 
The interested reader is referred to Chan dar et al.| q201 1) > 
for a description of the methods used to select clusters and 
to derive physical quantities such as their age and mass. 
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Fig. 7. — Properties of the GMC (left column) and young (< lOMyr) stellar cluster (right column) populations in different M5 1 environments, 
compared to the shape of the CO PDFs. The Spearman rank correlation is indicated at the top left of each panel. 



The relationship between CO emission and other SF tracers 
within the PAWS field is discussed in several companion pa- 
pers (Meidt et al., submitted, Schinnerer et al., in preparation). 

As for GMCs, we find evidence for a strong connection 
between the prevalence of bright CO emission and M51's 
young (t < 10 7 Myr) cluster population. In particular, more 
positive IDI values are associated with a higher maximum 
young cluster mass M VC 95 (defined analogusly to M gmc ^), 
and with a higher number surface density J\f yc and combined 
mass surface density M. yc of young clusters (Figure l7Te] 



to [g]). The origin of these trends would seem to lie in a phys- 
ical - as opposed to statistical - connection between young 
clusters and GMCs: M yc ,95, JV yc and M. yc are better correlated 
with the average GMC mass surface density ((Sh 2 )) and 
maximum GMC mass {M gmc ^) than with the total num- 
ber of young clusters or GMCs within each M5 1 environment. 

An exception to the good correspondence between the shape 
of the CO PDFs and the properties of GMCs and young 
stellar clusters is the slope of the cluster mass distribution: 
while there is a clear trend for the GMC mass spectrum to 
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become shallower in regions where bright CO emission is 
more prevalent, a connection between the slope of the young 
cluster mass function and the shape of the CO PDFs is less 
obvious (cf panels [d] and [h] of Figure [7] There is some 
indication that the mass distribution of the young cluster 
populations in the arm and interarm regions follow the same 
trend with IDI as GMCs, but within the central kiloparsec 
of M51 (i.e. the molecular ring and nuclear bar regions), 
the young cluster mass distributions are steep (< -2.5) even 
though bright CO emission is relatively dominant there. 

In Figure [8] we plot the slope of the young cluster mass 
spectrum directly against the slope of the GMC mass spec- 
trum for the different M51 environments. Various techniques 
for estimating the slope of the mass spectrum have been 
used by empirical studies of young clusters and GMCs 
(e.g. differential versus cumulative mass distributions, bins 
of equal width versus bins containing an equal number of 
objects), and the derived slope is known to be sensitive to 
factors such as cloud decomposition algorithm, the adopted 
low-mass completeness limit, and undersampling and/or the 
existence of a physical truncation to the distribution at high 
masses. We constructed both differential and cumulative 
mass distributions for the GMC and young cluster popula- 
tions, and for each variant, we estimated the slope multiple 
times using different binning strategies (in the case of the 
differential mass distributions) and mass ranges to calculate 
the fit. We describe these tests more fully in Appendix [C] 
The two panels of Figure [8] represent the results using a 
cumulative (panel [a]) and differential representation (panel 
[b]) for the mass distributions. In both panels, the error bars 
reflect the dispersion in the estimated slopes of the GMC and 
young cluster mass spectra in each environment. Despite the 
systematic uncertainties that limit the accuracy of any indi- 
vidual measurement of the mass distribution slope, we can 
therefore confidently draw two conclusions from Figure [8] 
The first is that while there is reasonable agreement between 
7 gmc and j yc for the arm and interarm environments of M5 1 , 
the slope of the mass distribution is not universal, i.e. the 
same values of 7 gmc and j yc do not hold everywhere within 
M51. Second, agreement between the slopes of the GMC and 
young cluster mass distributions is not ubiquitous. Averaged 
across the entire PAWS field and for some of the arm and 
interarm regions, 7 gmc and y yc agree to within ~ 0.3 dex, but 
for the molecular ring and upstream environments 7 gmc 56 7^, 
regardless of the method used to represent the mass functions. 

The trends in Figure [8] are remarkable since the observed 
similarity between j gmc sa -1.7 and 7 VC w -2.0 is frequently 
cited as evidence for the weak mass dependence of both 
the efficiency of star formation in GMCs and probability of 
cluster disruption. In their investigation of stellar feedback 
and disruption of GMCs, for example, |Fall et aT] (2010) 
derive relations between j gmc and j yc , which for isolated, 
bound systems in the absence of magnetic support are linked 
via the slo pe of the mass versus size relationship of GMCs. 
Fall et al. (2010) examine the limiting regimes of energy- and 
momentum-driven feedback, arguing that the constant surface 
density of GMCs (M ex R 2 ) ensures that 7 gmc ~ j yc regardless 
of the type of feedback that dominates GMC disruption. 

In a globally-averaged sense, the young cluster and GMC 
populations of M51 would seem to conform to the model 



outlined by Fall et al. (2010 1. In this case, 7 gmc = -1.7, 
7y C = -2. 1 and tor momentum-driven feedback, the predicted 
siope of the GMC mass-size relation is a = 1.8, in acceptable 
agreement with the observed value (a i, s = 2.1, Colombo 
et al., submitted). However, our results for individual 
environments within M51 suggest a more nuanced inter- 
play between molecular gas, young clusters and galactic 
structure. Downstream of the spiral arms, for example, 
Igmc ~ lye ~ -2.0 and a bs = 1 -8, in good agreement with the 
model prediction for energy-driven feedback. Yet upstream 
of the spiral arms (where 7 gmc ~ -f yc « -3.0), in the molecular 
ring (7 gmc ~ -1.4, 7,, c sa -2.7) and in the first spiral arm 
pattern of M51 (~y gm c ~ —1.5, j yc sa -2.1 ), the exponent of 
the mass-size relationship predicted by Fall et al. (2010 1 
is too shallow compared to the observed value by 0.4 to 
1.0 dex, regardless of the assumed feedback mechanism. 
This suggests that there may be regions within galactic disks 
where physical p rocesses that would seem to be excluded by 
Fall et al. (2010[), e.g. cluster coalescence, mass-dependent 



cluster/GMC disruption or a dominant role for energy-driven 
feedback, are in fact important. More generally, however, 
it highlights how valuable information can be lost in the 
calculation of galaxy-wide averages. With datasets that 
yield statistically significant samples of GMCs and other 
star-forming phenomena within ^kiloparsec-scale regions, 
it is timely that physical quantities (for example, based on 
a consideration of galaxy dynamics and/or ISM properties) 
determine the environments where GMC properties and 
extragalactic star formation are investigated, rather than 
relying solely on radial profiles and/or apertures that are 
'blind' to their location with respect to galactic structure. 



5.4. Comparison to Numerical Simulations of Galactic 

Disks 

Our analysis in this paper was prompted, in part, by numer- 
ical simulations showing that the gas density distribution 
in galactic disks is well-represented by a single l ognormal 
function spann ing several orders of magnitude (e.g. Wada & 
Norman 2007 ). If this is an accurate description of real galac- 
tic disks, then the result offers insight into the origin of the 
KS law, which can be reproduced from a lognormal density 
PDF with a limited number of plausible assumptions , such 
as a c ritical density threshold for star formation (Elmegreen 
2002), or a direct proportionality between the lo cal gas and 



star formation rate densities (e.g. |Kravtsov|20 03). While the 
overall /(CO) and T mb PDFs for the PAWS field are roughly 
lognormal, our results in Section |4.3| show that this average 
PDF shape obscures considerable diversity among the PDFs 
observed for different ~kiloparsec-sized regions within M5 1 . 
Since the regions that we use to investigate the PDFs are 
defined according to dynamical criteria, our basic result is 
that large-scale dynamical processes in M51's inner disk 
have an observable effect on the density (and column density) 
distribution of M51's molecular ISM. 

There are several further characteristics of our observed 
CO PDFs that are noteworth y in relation to the simulation 
results. For example, Wada & Norman (2007) find an 



increase in the logarithmic width of the PDF of ^0.3 dex 
for an order of magnitude increase in the mean gas density. 
This is roughly consistent with the variation in the width of 
the /(CO) PDFs of M51 and the LMC, between which the 
average H2 column density also varies by a factor of ~ 10 
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Fig. 8. — Slope of the young cluster mass function versus the GMC mass function for different dynamical environments within M51. Panel 
[a] represents the results using a cumulative formulation for the mass distribution; panel [b] shows the results for a differential representation. 
In both panels, the error bars reflect the dispersion in the estimated slopes (see AppendixO. The dashed diagonal line indicates equality, and 
the dotted diagonal line represents "/ gmc = y yc + 0.3 The plot annotations are the same as in Figure!/] 



(see Section |44| i. However, we caution that several effects 
limit the extent to which we can compare our observed CO 
PDFs to the density PDFs from simulations. The first is 
that the simulations describe gas density across six orders 
of magnitude, corresponding not only to the molecular ISM 
but also to the atomic and warm ionized phases. In practical 
terms, this makes the numerical result difficult to verify, since 
different observational tracers must be used to probe different 
phases of the interstellar gas, and each of them are sensitive 
to a much narrower range of densities (and column densities) 
than the full dynamic range of the simulated PDFs. 

The few simulations that include explicit treatment of molec- 
ular chemistry tend to find a range of densities between ~ 1 
and 10 3 cm" 3 for the H2 gas in galactic disks, and that the dis- 
tribution exhibits a sharp cut-off below the H2 self-shielding 
limit (n < 5 cm" 3 , see e.g. figure 11 of |Dobbs et al.||2008| l. 
The use of 12 CO(7 = 1 — >• 0) emission to trace the H2 column 
density should narrow the observed PDF even further, since 
at moderately low extinction (Ay ~ 1 mag, e.g. |Wolfire et al. 
2010) H2 can self-shield while CO molecules are photodisso- 
ciated. At high H2 column densities (N(H2) > 10 22 cm" 2 ), on 
the other hand, the CO-emitting structures within a GMC will 
start to overlap and shadow eac h other, lea ding to a saturation 
of /(CO) in tensities (e.g. |Shetty et ap[201l) . |Feldmann| 
|et al.| ( 2012) 1 show that this saturation should occur at /(CO) 
intensities near a few 100 Kkms" 1 . In general, however, 
the /(CO) PDFs in Figure HI show no evidence for a peak at 
these intensities caused by a 'pile-up' of saturated pixels, 
but instead more closely resemble their pure N{H2) PDFs 
after they exclude pixels with l ow CO integrated inte nsities 
(< 0.2 Kkms" 1 , see figure 8 of |Feldmann et al. |20 1 2| . We 
suggest that this is because the typical CO linewidths in M51 
are much greater than the two possibilities considered by 



Feldmann et al. (2012 1 (i.e. a constant linewidth of 3 km s 



and hence /(CO) remains a relatively good tracer of the Fb 
column density. This is consistent with recent studies of gas 
and dust in M51, which suggest that the Xqo f actor has a 
roughl y Galactic value th roughout the disk (e.g. |Schinnerer| 
|et al.|2010||Tan et al.|20TT]|Mentuch Cooper et al.|2012) . 

What insights can models provide regarding the diversity 
of distribution shapes that we observe? As noted in the 
Introduction, the PDF for supersonically turbulent isothermal 
gas is lognormal when the influence of gravity is negligible. 
In this case, the logarithmic width of the PDF x varies with 
the Mach number M . according to x 2 ca ln(l + 0.25.M 2 ) 
(Pado an~et al.| [1997). Across a galactic disk, however, the 
temperature and average density of the molecular gas will 
vary with location. Thus, a more realistic expectation for 
the observed density PDF on global to kiloparsec scales may 
be the convolution of the local lognormal PDF (reflecting 
the distribution of densities within a region over which the 
average gas density p ave and Mach number are relatively con- 
stant) with the PDFs of /9 ave and Ai within the galactic disk. 
Elmegreen (|2011) has recently considered such a model, 



presenting convolution PDFs for several idealized cases of 
gas clouds with different radial density profiles, and variable 
Mach numbers (see his figure 1). The resulting PDFs clearly 
diverge from a pure lognormal shape. As clouds become 
more centrally condensed (and hence more dominated by 
self-gravity), the convolution PDFs develop a power-law 
tail at high densities, with a slope that varies inversely 
with the slope of the density profile. If the Mach number 
decreases at higher average densities, on the other hand, the 
convolution PDF appears truncated relative to a pure log- 
normal since the local PDFs get narrower with increasing p ave . 



or virial scaling of the linewidth with mass surface density). 
As a consequence, the majority of the CO-emitting molecular 
structures in M5 1 do not shadow each other in velocity space 



In broad terms, the analysis by Elmegreen (2011) suggests 
that the diverse shapes of the /(CO) PDFs in Figure fflreflect 
large-scale variations in the average density, temperature 
and/or velocity fluctuations for the molecular gas within 
different M51 environments. In reality, these properties are 
likely to vary simultaneously, so attributing a specific PDF 
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morphology to a variation in one physical quantity and/or 
process is not straightforward. Nevertheless, it is suggestive 
that the /(CO) PDFs in the interarm region - where we might 
expect the temperature, density and velocity structure of the 
molecular gas to be determined by cloud-scale processes 
- resemble the pure lognormals expected for isothermal 
supersonically turbulent molecular gas, whereas the /(CO) 
PDFs in M51's spiral arms - where the molecular gas not 
only reaches higher densities, but its velocity structure 
can be influenced by large-scale dynamical effects such as 
streaming motions and the spiral shock - more obviously 
diverge from a lognormal shape. The molecular ring region, 
where the /(CO) PDF is very broad and almost flat-topped, 
is arguably the extreme case both in terms of dynamical 
effects and star formation activity. Although shear and 
large-scale non-circular motions should be low in the ring, 
the molecular gas accumulates here due to opposing bar 
and spiral torques and the average gas velocity dispersion is 
relatively high (Colombo et al., submitted). The level of star 
formation activity in the ring is also high (Schinnerer et al., in 
preparation), so feedback from nascent stars may also have a 
strong effect on the distribution of gas densities in this region. 

Finally, we note that Hopki ns" et al.| ((2012) have recently 
shown that the dominant mode of stellar feedback (and not 
just total amount of star formation) has an observable effect 
on the shape of the density PDF for the cold gas component 
in their simulated galaxies. Like Wada & Norman (2007), 
they find that the width of the density PDF decreases for sys- 
tems with lower average gas densities. The density PDFs 
of their simulated gas disks show striking departures from 
lognormality (see their figures 10 and 11), however, which 
they attribute to the inclusion of cooling, self-gravity and 
a physically-motivated implementation of different feedback 
mechanisms in their simulations. In particular, they find that 
radiation pressure is crucial for suppressing a pile-up of gas 
with high densities (n > 10 4 cm" 3 ), since pure gas heating 
(e.g. by supernovae, stellar winds and Hn photoionization) 
is ineffective in disrupting dense gas clumps where the cool- 
ing time is much shorter than the dynamical time. While the 
overall shapes of the /(CO) and T mh PDFs for different M51 
environments almost certainly reflect the combined action of 
several distinct physical processes, the absence of a secondary 
peak in the PDFs at high CO intensities would seem to con- 
firm that the dominant feedback mechanism in M5 1 must be 
effective at preventing the build-up of high-density material. 

6. CONCLUSIONS 

In this paper, we have presented the probability distribution 
functions (PDFs) of CO integrated intensity and CO bright- 
ness within the inner disk of M51, using new high resolution 
(~ 40 pc) data from the PdBI Arcsecond Whirlpool Survey 
(PAWS, Schinnerer et al., in preparation). We have compared 
the PDFs of these properties for different environments 
within the PAWS field, and to PDFs constructed using high 
resolution CO datasets for two nearby dwarf galaxies. We 
report the following results and conclusions: 

1. On ~ 40 pc scales, the distribution of /(CO) integrated 
intensities within the inner ~ 1 1 x 7 kpc of M5 1 spans ~ 1 .5 
orders of magnitude above the 3a sensitivity limit of the 
PAWS data. The shape of the /(CO) PDF is consistent with 
a lognormal function with a mean of 20 Kkms" 1 and a 



logarithmic width of 0.4. Relative to this LN function, there 
is some evidence that the observed PDF is truncated for 
/(CO) values greater than ~ 200 Kkms -1 . 

2. The CO brightness temperatures that we measure for the 
inner disk of M51 span ~ 1 to 10K, where the lower limit 
corresponds to our survey's 3a sensitivity limit. The shape of 
the T m \, PDF can be represented by a lognormal function with 
a mean of 0.9 K and a logarithmic width of 0.3, but a broken 
power law with a slope of ~ —1.4 for 1 < T m \, < 5K and a 
much steeper slope of ~ -4.9 for T m ^ > 5 K is an equally 
adequate description of the distribution. 

3. The CO PDFs that describe the emission in the inner 
disk of M51 are clearly different to the PDFs obtained for 
M33 and the Large Magellanic Cloud (LMC). The maximum 
/(CO) and T m b values observed in M5 1 are 1 to 1 .5 dex higher 
than in the two low-mass galaxies. The CO PDFs in M51 are 
also wider, consistent with numerical results indicating that 
the width of the density and column density PDF in creases 
with the avera ge gas density of a galactic disk (e.g. Wada & 
|Norman|2007) . 

4. The CO PDFs for different dynamical environments within 
M51's inner disk exhibit diverse shapes. The CO PDFs in the 
interarm regions are narrower than in the spiral arms, nuclear 
bar and molecular ring regions. The distributions of /(CO) 
and r m b are approximately lognormal in the interarm, while 
the PDFs in the arms, ring and bar exhibit strong departures 
from lognormality such as power-law slopes and/or trun- 
cations at high CO intensities. While a lognormal function 
may provide an adequate description for the overall gas 
distribution within a galaxy, phenomena such as streaming 
motions, spiral arm shocks and star formation feedback 
produce observable changes to the gas density distribution 
for ^kiloparsec-sized within galaxies. 

5. To avoid assuming a particular functional form for the 
CO PDFs in M51, we characterised their shape using the 
brightness (or integrated intensity) distribution index, origi- 
nally devised by |Sawada et al.| ( |20l 2b), a simple parameter 
that specifies the ratio between Dright and faint emission. 
With this, we showed that the shape of the CO PDFs for 
dynamically-defined, kiloparsec-scale environments within 
M51 are strongly correlated with physical properties of 
the GMC and young stellar cluster populations of those 
environments and, we infer, their star formation activity. The 
implications of this result for interpreting the observational 
scatter in extragalactic star formation laws are explored in 
several companion papers (Meidt et al. submitted, Schinnerer 
et al., submitted, Leroy et al., in preparation). 

6. Consistent with the predictions from numerical simulations 
(e.g. Wada & Norman 2007), we find a shallow increase in 
the width of the PDF with increasing average gas surface 
density. The dynamic range of the observed /(CO) PDFs is 
also in approximate agreement with the distributions of H2 
column density obtained by simulations that include explicit 
treatment of molecular chemistry (e.g. |Dobbs et al.||2008] 
Feldmann et al. 2012), but we do not observe a secondary 
peak in the PDFs at high CO intensities corresponding to CO 
saturation. We suggest that this is because the CO linewidths 
in M5 1 are typically larger than the linewidths adopted by the 
simulations, so the 'mist-model' explanation (Dickma n et al.| 
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1986) of the Galactic Xco factor remains valid even at high 



H2 column densities. 

7. We show that the diverse shapes of the CO PDFs in M51 
are qualitatively similar to the deviations from lognormal- 
ity expected from the combined action of star formation 
feedback and large-scale variations in density, temperature 
and velocity structure throughout M51's inner disk. Our 
results suggest that star formation feedback on small scales 
and dynamical effects on large scales (e.g. the influence 
of the stellar bar and spiral density way) together regulate 
the velocity structure of the molecular gas, and that these 
processes in combination with gas self-gravity determine 
the shape of the CO PDFs. Isolating the dominant physical 
process responsible for the morphology of each PDF will 
require a more detailed comparative analysis with theoretical 
models however. 

8. The precise shape of the /(CO) and T m t, PDFs is sensitive to 
several non-physical effects including resolution, sensitivity, 
and the method used to identify significant emission within a 
spectral line cube. These caveats should be kept in mind by 
future studies that compare the PDFs derived from CO ob- 
servations of different galaxies, or aim to validate numerical 
models using observational results. In particular, we note that 
the estimated logarithmic width of the PDFs tends to decrease 
for datasets with poorer sensitivity, and that degrading the 
observational resolution to a spatial scale greater than the 
characteristic spacing between high brightness structures can 
produce the appearance of a threshold in the PDFs. This 
result is not unique to PDFs of CO emission and suggests 
that observations of thresholds on kiloparsec scales should be 
interpreted with care. 
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APPENDIX 
MASKING METHODS 

As discussed by Pety et al. (submit ted), a number of different techniques for constructing CO integrated intensity images have 
been presented in the literature (e.g. Heifer et al.|2003 Dame|20TT Wong et al.|201 1 1. These include: 



1. A sigma- clipping method, Ml, whereby pixels containing emission below tiurms are blanked, grms is the RMS of the 
noise variations, which we calculate for each independent line-of-sight. For the comparison in this Appendix, we adopted 
re = 3. 

2. A dilated mask method, M2, which identifies islands of significant emission by selecting peaks above a threshold of togMS 
across two contiguous velocity channels. The preliminary mask is then expanded to include all contiguous pixels with 
emission above eajwis- We adopted (t,e) = (5, 1.2). 

3. A smooth-and-mask method, M3, which generates a version of the cube that has been spatially smoothed to an angular 
resolution of 6, and identifies emission in the smoothed cube above a significance threshold iikjrms- The blanking mask is 
then transferred back to the original (i.e. full resolution) data cube. We adopted {0,m) = (3"6,5). 

4. An Hi velocity prior method, M4, which assumes all of a galaxy's CO emission arises in velocity channels within a 
restricted interval, AV , around the radial velocity corresponding to the peak of the H i line profile for each line-of-sight. 
We used AV = 50 km s" 1 . 

In addition to these, we defined a final mask for the PAWS cube that optimized flux recovery while eliminating anomalous 
features in the map of M51's velocity field (M5). The construction of this mask is described in detail by Colombo et al. (in 
preparation). Example /(CO) maps for the PAWS field constructed using each of the five techniques are shown in figure 23 of 
Pety et al. (submitted). 

In Figure 151 a], we show the PDFs obtained from the different /(CO) maps. It is clear that the map construction method affects 
the shape of the /(CO) PDF. Differences between the curves are apparent up to ~ 60 K km s" 1 , which is considerably greater than 
our nominal 5a RMS sensitivity limit (^18 Kkms" 1 ). The mean (i ) and logarithmic width (x) of the best-fitting LN function to 
each PDF in Figure |5Ja] are listed in Table [3] The PDF from the /(CO) map constructed using the dilated mask (M2) is the most 
similar to the PDF from our preferred mask (M5), although it recovers fewer pixels than M5 with 10 < /(CO) < 60 Kkms -1 . 
The same applies to the PDF corresponding to the H i velocity prior method (M4), although this PDF also slightly underestimates 
the number of high intensity pixels (/(CO) > 160 Kkms -1 ). Inspection of the CO velocity dispersion map for the PAWS field 
shows that there is a small fraction of pixels with FWHM linewidths greater than 17 kms" 1 (mostly in the nuclear bar and 
southern spiral arm regions), so our chosen velocity interval of 50 km s" 1 excludes some genuine emission in the wings of these 
line profiles. The smooth-and-mask technique (M3) recovers the least emission at intermediate /(CO) values. This is because 
we used a relatively large smoothing kernel and high significance threshold, so compact regions with moderate significance in 
the original data cube are excluded from the final map. The large smoothing kernel is also the reason why the M3 map does 
not show a sharp cut-off at a low /(CO) value, as it incorporates many pixels with low significance that are adjacent to high 
brightness regions. Unlike the other methods, which tend to peak around /(CO) ~ 20 kms" 1 , the sigma-clipping method (Ml) 
peaks at our nominal 3ctrms sensitivity limit. Even though we define ctrms locally, i.e. we estimate vrms for each line-of-sight, 
inspection of the /(CO) map constructed using Ml reveals that many of these pixels come from the edge of the field where the 
brightness sensitivity of PAWS declines. Although this masking method is relatively common, we regard the resulting PDF to be 
the least reliable measure of the /(CO) distribution in M5 1 's inner disk. In summary, while the different methods for constructing 
the /(CO) map yield PDFs that are different in detail, all the /(CO) PDFs except that obtained using Ml are roughly consistent 
with a LN function with mean so ~ 20 Kkms" 1 and logarithmic width x = 0.4, which we have indicated by a thick cyan line 
in Figure I5J a] (see also Table [3]). Future studies should keep in mind that estimates for the shape of the /(CO) PDF are more 
likely to be dominated by systematic uncertainties due to different techniques for identifying significant emission than by simple 
counting and/or measurement errors. 

As well as the /(CO) PDF, we tested the effect of different masking techniques on the shape of the T m ^ PDF. The results for 
the central 140" x 90" of the PAWS field are presented in Figure [9jb]; we do not construct PDFs using the entire field since the 
noise increases significantly towards the edge of the map (see Figure [T]) which makes the resulting PDFs harder to interpret. For 
completeness, we also show the T m \, PDF of the PAWS cube without applying any mask (M0, grey histogram). Even though the 
total CO flux of the different masked and unmasked cubes agrees to within ~ 30% (see table 8 of Pety et al., submitted), the grey 
and magenta histograms (M0 and Ml) only converge with the other PDFs for r m b > 2 K. Below 2 K, the M0 and Ml histograms 
begin to curve upwards, departing from the roughly LN shape of the distribution at higher intensities. The remaining PDFs are 
more similar, suggesting that the masking techniques that use additional criteria (e.g. proximity to a bright peak) are successful 
at retaining genuine emission at ~ 3<trms, whereas a simple 3cr clip (i.e. masking method Ml) may retain a significant number 
of isolated noise peaks. Never the less, some genuine low brightness emiss ion co uld b e masked by methods M2 to M5, although 
its spatial distribution is difficult to determine. In our analysis (Sections |4. 1 1 to [4~4~]>, we focus on the shape of the /(CO) and 
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Fig. 9. — PDFs of [a] /(CO) and [b] r m b for the PAWS field. The different colours represent PDFs obtained from different masking techniques 
to identify significant emission within the data cube (see text). A dashed parabola indicates the LN function that provides the best-fit to the 
corresponding PDF. The thick cyan line in panel [a] represents a LN function with mean (/(CO)} = 20 Kkm s and logarithmic width x = 0.45; 
in panel [b], the cyan line represents a LN function with mean (r m b) = 1.3 K and logarithmic width x = 0.35. In both panels, the grey shaded 
region represents values beneath the nominal 3/jrms sensitivity limit. The top horizontal axis indicates the equivalent PL surface density for the 
/(CO) or T ra b value on the lower axis, assuming Xco = 2 x 10 20 cnT 2 (Kkms~')~' and a helium contribution of 1.36 by mass. The error bars 
represent the uncertainty associated with simple counting (VN) errors. 



TABLE 3 
Fit Parameters for CO PDFs in Figures[91and[T21 



Figure CO Property Description 



Mean 
•So 



Logarithmic Width Goodness-of-fit 



9 a] 


/(CO) 


9 a] 


/(CO) 


9 a] 


/(CO) 


9 a] 


/(CO) 


9 b] 


Zmb 
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/mb 


9 b] 


/inb 
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12 a] 


/(CO) 


12 a] 


/(CO) 


12 a] 


/(CO) 


12 a] 


/(CO) 
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/(CO) 
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/(CO) 


12 c] 


/(CO) 



Ml 

M2 

M3 

M4 

MO 

Ml 

M2 

M3 

M4 

<r ms = 0AK,(t,e) = (5,1.2) 

cr RM s = 0.6K,(t,e) = (5,l.2) 

<r ms = 1.0K,(t,e) = (5,1.2) 

a ms = 2.0K,(t,e) = (5,1.2) 

<r«s = 0.4K,((,e) = (3.5,2) 

<r ms = 0.6K,(t,e) = (3.5,2) 

(TMS = 1.0K,(f/) = (3.5,2) 



6.3 Kkms-' 


0.62 


0.10 


15.3 Kkms- 


0.50 


0.08 


21.8 Kims" 


0.45 


0.07 


24.4 Kkm s" 


0.36 


0.05 


0.9 K 


0.40 


0.19 


0.9 K 


0.40 


0.18 


1.2K 


0.36 


0.20 


I.4K 


0.34 


0.19 


I.4K 


0.34 


0.20 


15.4 Kkms~ 


0.50 


0.18 


21.0 Kkms~ 


0.44 


0.09 


22.1 Kkms- 


0.43 


0.14 


62.0 Kkm s" 


0.31 


0.14 


10.1 Kkms~ 


0.66 


0.13 


21.1 Kkms~ 


0.39 


0.12 


28.9 Kkms" 


0.33 


0.11 



Parameters of best-fitting functions to PDFs in Figure R)] The parameters of the LN functions are determined 
from a Levenberg-Marquardt fit to Equation[3] We use tne logarithmic dispersion of the fit residuals to estimate 
the goodness-of-fit. 



r m b PDFs at relatively high brightness (i.e. brighter than 4<7[>ms), so our results and interpretation should not be affected by the 
presence of such a faint emission component. 



As for the /(CO) PDF, the smooth-and-mask method (M3) recovers the least emission at intermediate r m b values (1 < 7" m b < 3 K) 
because compact, high brightness regions are diluted beneath our 5vrms threshold in the smoothed cube. The dilated mask 
technique (M2) yields a PDF that appears more like a broken power-law with a turnover at ~ 3 K. This variation is reasonable, 
since M2 should exclude some isolated regions of genuine emission that fall beneath our 5<trms threshold, but include some 
noise at the edges of the mask with T m \, ~ 2<7rms- For values above ~ 3<trms, the T m \, PDFs for M2, M4 and M5 are practically 
identical, exhibiting a mean T m \, ~ 1 .3 K and logarithmic width x ~ 0.35. 
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Fig. 10. — The PDFs of [a] /(CO) and [b] T m b within the PAWS field, obtained after convolving the original PAWS data cube with Gaussian 
smoothing kernels of varying width (see text). The PDFs are obtained from cubes where significant emission is identified using a dilated mask 
method. The error bars represent the uncertainty associated with simple counting (yN) errors. 



We use our tailored m ask, M5, for our analysis of the CO emission in the PAWS field (Section 4. 1 1 and within different M5 1 
environments (Section pO) . As it closely reproduces the M5 results for M51, but is much simpler to implement acr oss multiple 
data sets, we use the dilated mask technique (M2) for our comparative analysis of M51, the LMC and M33 (Section |4~4| ). 

RESOLUTION AND SENSITIVITY EFFECTS 

Previous studies of the column density PDF for individual molecular clouds have shown that the shape of the PDF depends 
on the spatial resolution of the data (e.g. Froebrich & Rowles 2010). To assess the importance of this effect, we smoothed the 
original M51 data cube to angular resolutions of 1.5, 3, 6, 12 and 24", corresponding to linear scales of 60, 1 10, 230, 450 and 
910pc respectively for our assumed distance to M51. We constructed /(CO) maps from all the cubes after applying the dilated 
mask method to identify significant emission. The resulting /(CO) PDFs are presented in Figure 10 a], while the PDFs of CO 
brightness are shown in Figure 10 b]. For the /(CO) PDFs, moderate variations in the resolution (Te. up to ~ 200 pc) produce 
PDFs with a similar shape over a restricted range of /(CO) values (20 to 100 Kkms" 1 ). It is striking, however, that the shape 
of the PDF at high intensities steepens as the angular resolution of the datacube is degraded. At ~ 1 kpc resolution, the /(CO) 
PDF appears to show a threshold at /(CO) ~ 60 Kkms" 1 ; smoothing the data over even larger scales has no further effect. A 
kiloparsec corresponds roughly to the radius of the central region, and also to half the distance between the spiral arms. The 
appearance of an upper threshold would therefore seem to occur because the emission from these high brightness regions is 
averaged together once the resolution is coarser than this spatial scale. 

We obtain similar, though not identical, results for the T m <o PDFs. The slope of the PDFs at intermediate intensities (i.e. from 
~ 5<Jrms to the r m b value where the PDF begins to steepen) appears relatively constant, flattening slightly as the resolution 
degrades from a slope of ~ -1.4 at 60 pc resolution, to ~ -0.7 at 450 pc resolution. We note that this is the opposite trend to 
what is observed for the /(CO) PDFs, which become steeper at lower resolution. Like the /(CO) PDFs, however, the T m ^ PDFs 
show a truncation that shifts to lower CO intensities as the smoothing scale increases. By ~ 1 kpc resolution, the PDF has a sharp 
cut-off at r m b ~ 1 K f or the same reason that the /(CO) PDF shows a threshold at this scale. The results of our resolution tests 
suggest that the appearance of thresholds and/or truncations in PDFs of gas emission tracers should therefore be interpreted with 
some caution, although these effects would appear to be most severe once the resolution of a dataset becomes comparable to the 
characteristic spacing between regions of high brightness (e.g. the spacing between spiral arms and/or between star formation 
complexes). Since the PAWS observations resolve these spatial scales, we do not regard resolution effects to be the main driver 
of the deviations from lognormality that we observe for the CO PDFs in Figures [4] and E\ The fact that the CO PDFs of M51's 
interarm regions do not show truncations at high intensities would further tend to support our interpretation that the truncations 
observed for the PDFs of the spiral arm and central environments are not solely due to limited spatial resolution. Since Galactic 
GMCs exhibit variations in their CO surface brightness, we would expect observations that spatially resolve internal structure of 
extrgalactic clouds to reveal more features in the shape of the /(CO) and T m \, PD Fs at high CO intensities (such as power-law 
tails due to the formation of strongly self-gravitating clumps, see e.g. Kainulainen et al. 2009). While resolution effects do not 
seem to be the primary explanation for the variation in PDF shap es with galactic environment that we describe in Section HI we 
therefore echo the recommendation by Wada & Norman (2007 see their figure 7) that well-matched resolution is critical for 
comparative studies between observational datasets, or between models and observations. 
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Fig. 11. — The PDFs of [a] /(CO) and [b] T ra b within the PAWS field, obtained after folding the original PAWS data cube along the spectral 
axis to have velocity channels of varying width. The PDFs are obtained from the resulting cubes after identifying significant emission using a 
dilated mask method. The error bars represent the uncertainty associated with simple counting (s/~N) errors. 



In addition to the effect of spatial resolution, we checked whether variations in the width of the velocity channels, i.e. the cube's 
spectral resolution, produced systematic changes in the shape of the /(CO) and T m \, PDFs. We constructed PDFs from cubes that 
had been folded along the velocity axis to 10, 20, 30 and 50 km s" 1 , and show the results in Figure 1 1 As the width of a velocity 
element increases, the /(CO) PDF narrows but the shape of distribution at high intensities (/(CO) > 50 Kkms" 1 ) remains 
unchanged. The T m \> PDF broadens with increasing channel width, but its shape is also relatively robust. The maximum observed 
CO brightness declines by ~ 3K as the channel width increases from 5 to 50 kms" 1 ; this is due to the spectral equivalent of 
beam dilution. 

In principle, the shape of the /(CO) and T m ^, PDF at high intensities should be relatively robust to noise, provided that the 
signal-to-noise ratio is sufficiently high. In practice, however, the sensitivity of most extragalactic 12 CO(7 = 1 — > 0) mapping 
surveys is limited. In Figure 12 a] and [b] we show the PDFs of integrated intensity and CO brightness, after adding increasing 
levels of Gaussian noise at the beam scale and using a dilated mask with (t,e) = (5,1.2) to identify regions of significant 
emission. For a moderate decrease in sensitivity {crms S [0-4, 1.0] K), the /(CO) PDFs retain their shape at high intensities (i.e. 
above ~ 50 Kkms" 1 ). Once the noise level is increased to <trms = 2K, however, the PDF begins to diverge significantly from 
the distribution obtained for the original cube, even at high intensities. This is because several regions in the nuclear bar and 
molecular ring regions are excluded from the mask once this noise level is added to the cube. High /(CO) measurements in these 
regions are often due to line profiles that are unusually wide (FWHM greater than ~ 30 kms" 1 ) but have moderate brightness 
in individual channels, so they do not possess a bright core that lies above 5<Jrms = 10 K. Even for moderate levels of added 
noise {a RMS < 1 K), however, the best-fitting LN functions to the /(CO) PDFs in Figure 12 a] and [c] appear to narrow as the 
sensitivity decreases (see Table 13). This occurs even though we restrict the range of CO intensities used to estimate the fit to 
pixels where the emission is brighter than ^urms- The dependence of the /(CO) PDF width on brightness sensitivity should be 
kept in mind when comparing the M51 PDFs to results from other galaxies or numerical simulations. 

The r m b PDFs from cubes with differing noise levels show more variation. Not surprisingly, the PDFs become narrower 
as the noise increases (since less emission satisfies our criteria for significance) but they also become steeper at low CO 
intensities, losing the appearance of a truncation at 7; n b ~ 5 K that is observed for the original cube. Since the construction of the 
dilated mask depends on the signal-to-noise, we examined whether similar trends were observed when we used different (t,e) 
combinations to identify significant emission. As an example, /(CO) and T m \, PDFs for varying levels of noise using a dilated 
mask with (t,e) = (3.5,2) are shown in Figure [T2|c] and [d] respectively. While these PDFs are not identical to the PDFs in 
panels [a] and [b], the effects of increasing noise on the shape of PDFs that we have described are not sensitive to the particular 
combination of (t,e) that we adopt to mask the input cubes. 



ESTIMATING THE SLOPE OF THE GMC AND YOUNG CLUSTER MASS FUNCTIONS 

For both GMCs and stellar clusters, the shape of the mass distribution is an important empirical signature of the physical pro- 
cesses that regulate their formation and disruption. Defined as the number of objects per unit mass, f(M) = dN/dM, numerous 
studies have f ound that the mass distributio n can be well-d escribed by a power law, f(M) ex M 13 , with a typical exponent of 
/3 s=i — 1.7 (e.g. Rosolowsky|2005 Fukui & Kawamura|2010 1 for extragalactic GMCs identified in 12 CO(7 = 1 — ► 0) surveys, and 
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Fig. 12.— PDFs of CO(7 = 1 -> 0) integrated intensity (panels [a] & [c]) and CO brightness (panels [b] & [d]) within the PAWS field, 
obtained after adding Gaussian noise to the datacube at the beam scale. The PDFs were obtained from cubes where significant emission was 
identified using the dilated mask method with (t,e) = (5, 1.2) (panels [a] & [b]) and (t,e) = (3.5,2) (panels [c] & [d]). In panels [a] and [c], a 
dashed parabola indicates the LN function that provides the best-fit to the corresponding PDF. The vertical dotted lines indicate our nominal 
3ct R ms sensitivity limits. The PDFs in the panels [b] and [d] were obtained using the central quarter of the PAWS field only, to suppress effects 
that are caused by the lower signal-to-noise at the edge of the PAWS field. The error bars represent the uncertainty associated with simple 
counting (yN) errors. 



(3 ss -2.0 for young star clusters (e.g. Chandar et al. 2010). An accurate determination of the shape of the mass distribution is 
crucial if it is to be used as a metric to quantify differences between and/or among populations of GMCs and clusters, and hence 
to argue for (or against) the universality of the processes that determine their evolution. As noted by several authors, however, 
it can be difficult to measure the shape of the mass distribution robustly, especially when the sample of observed objects is 
small. Well-recognized sources of uncertainty include the method used to identify structures of interest (which usually depend 
on both resolution and sensitivity), unambiguous determination of the low-mass completeness limit, inadequate sampling of the 
high-mass end of the distribution, uncertainty in mass measurements of individual objects, and (for mass distributions that model 
a differential formulation with a histogram) the choice of binning parameters. An added complication is that there are several 
'standard' methods for representing mass distributions employed by the GMC and stellar cluster research communities. In broad 
terms, studies of stellar cluster populations tend to adopt a differential formulation of the mass distribution, either separating 
the cluster mass measurements into bins of variable width with equal numbers of clusters in each bin, or into bins of uniform 
logarithmic width but with a variable number of clusters in each bin, while it has become increasingly common for studies of 



CO PDFs in M51 by PAWS 



27 






5 


6 


7 




1 


i 


Ml 

I : 

ring 

iS 

I I 


i 1 i 
i 

i 
i 
i 




1 


i 


1 1 

!' 

i! 

i i il ii 


IB, 

! 
i 
1 1 1 







s = log(Mass/[M so i]) 



Fig . 13. — Cumulative mass distributions for GMCs in the different M5 1 environments. In each panel, the vertical dot-dot-dashed lines indicate 
the mass range that was used to obtain the fit in the example shown. The black dashed lines indicate the limiting mass range over which we 
conduct the trials. 



GMC populations to represent the GMC mass distribution in cumulative form. 

In light of these uncertainties and variations in technique, it is evident that any literature comparison between the parameters 
of mass distributions must be made with caution. In this Appendix, we estimate the uncertainty in the slope of the GMC and 
stellar cluster mass functions for different M51 environments by constructing the distributions using three common methods, 
and by using different mass ranges in the linear regression that we use to fit the mass distribution. Our analysis falls short of 
being a general study of bias in mass distributions in at least two important ways. Firstly, we model all the mass distributions as 
pure power-laws, and do not investigate other functional forms (such as a truncated power-law or Schechter function) that have 
been used to describe GMC and cluster populations. Secondly, we only examine a simple estimator - i.e. ordinary least squares 
linear regression - to deter mine the best-fitting power-law to the mass functions. Since our main goal is to assess whether the 
trends discussed in Section 5.3 are robust, we consider it sufficient that we have used the same power-law model a nd stat istical 
estimator f or all the G MC and cluster mass distributions in our analysis. We refer the reader to other studies (e.g. Rosolowsky 
2005 1 Reid et al .|2010[ ) for more comprehensive investigations of the general problem of estimating the true mass distribution of 
objects from observational data. 



In Figures 13 to 16 we plot example mass distributions for the GMCs and young (< lOMyr) clusters in the eight M51 environ- 



ments that we examine in this paper. The distributions in Figures 13 and 15 are constructed using a cumulative representation, 
while those in Figures 14 and 16 are differential mass distributions for the same GMC/cluster populations, constructed using one 
example set of binning parameters. In each panel of Figures [14] and [16] we show both common forms of the differential mass 
distribution, i.e. a histogram with bins of equal logarithmic width (grey lines) and a histogram with equal number of objects 
per bin (black points). The black vertical dashed lines indicate the limiting mass range over which we estimate the fit for each 
distribution; the limits of this range were chosen so that the fit was calculated over an appreciable range of objects masses but 
avoiding regions of the mass distributions that show clear evidence for incompleteness effects (i.e. flattening) at low masses or 
truncation and/or sampling effects at high masses. The grey dot-dot-dashed lines indicate the actual mass range that was used to 
obtain the fit in the examples shown. 

The binning parameters and mass limits that were used in the linear regression for our trial fits to the mass distributions 
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are listed in Table |4| For most environments, we tested mass ranges with lower (upper) limits logM, 
(logM, 



G [5.8,6.1] 
€ [6.5 ,"0]) for GMCs and logM, m „. e/ e [3.6,3.9] (logAf mflA -. e / e [4.2,4.4]) for young clusters. Additionally, for the 



min,gmc 



stellar clusters, we required that the mass range used to estimate the fit was larger than 0.4 dex, i.e. logM mflA > logM„„„ + 0.4. 
For the fits to the differential mass spectra, bins containing less than two objects were excluded from the fit, and the fit 
was only estimated when three or more bins occupied the specified mass range. For the ring environment, the range of 
lower mass limits for stellar clusters was modified to logM m ,„ c / e [3.75,3.9] because the mass distribution (see panel [b] 
of Figure 15 1 flattens sharply below M ~ 6000 M . In principle, this flattening could have an observational (e.g. a higher 
completeness limit due to crowding/extinction in this region) or a physical origin (e.g. a higher probability of disruption for 
low-mass clusters objects). These possibilities will be explored in a future paper that investigates the shapes of the GMC 
and cluster mass distributions in detail; here, we simply reduce the mass range that we fit to the part of the distribution that 
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Fig. 14. — Differential mass distributions for GMCs in the different M51 environments. In each panel, the grey lines indicate a histogram with 
bins of equal logarithmic width, while the black points represent a histogram with equal number of GMCs per bin. Other plot annotations are 
the same as in Figure[T3] 
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Fig. 15. — Cumulative mass distributions for young (< lOMyr) stellar clusters in the different M51 environments. Plot annotations are the 
same as in Figure[T3] 



conforms to a pure power-law. In total, we obtained 35 (49) estimates of the cluster (GMC) mass distribution slope for each 
of the cumulative mass distributions, and between 36 (362) and 385 (539) estimates for each of the differential mass distributions. 



18 For both 



Histograms showing the distribution of slopes that we obtain from the trial fits are shown in Figures 17 and 
GMCs and young clusters, it is evident that the slopes obtained using a cumulative representation of the masscTistributions 
are systematically steeper (more negative) than for the differential formulation, although the offset between the peak of the 
histograms varies with environment. The discrepancy likely reflects the fact that the observed mass distributions are not pure 
power-laws, but tend to steepen continuously across the observed range of masses. This steepening is not well-captured by the 
differential distributions, since the values in each bin are weighted towards the lower-mass objects (which are more common) 
and hence a systematically shallower slope. The agreement between the two techniques is better for environments (e.g. upstream 
of the spiral arms) where the mass distribution more closely follows pure power-law behaviour, consistent with this interpretation. 



CO PDFs in M51 by PAWS 
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Fig. 16.- 
Figure[T4"| 
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s = log(Mass/[M so i]) 

Differential mass distributions for young stellar clusters in the different M51 environments. Plot annotations are the same as in 



TABLE 4 
Parameter Space Explored By GMC/Cluster Mass Distribution Trials 



Parameter 



Range 



Comment 



n Un [8,9, 10, 11, 12, 13, 14, 15, 16, 17, 18] 

n obj [8,9, 10, 11, 12, 13, 14, 15, 16, 17, 18] 

logM^rf" [3.6, 3.65, 3.7, 3.75, 3.8, 3.85, 3.9] 

logM mn ,. c , [4.2,4.25,4.3,4.35,4.4] 

logM m ,„. gmc [5.8,5.85,5.9,5.95,6.0,6.05,6.10] 

\ogM max , gmc [6.5,6.55,6.6,6.65,6.7,6.75,6.8] 



Number of histogram bins 

Number of GMCs/cluster in each bin 

Lower mass limit for fit to cluster mass distribution 

Upper mass limit for fit to cluster mass distribution 

Lower mass limit for fit to GMC mass distribution 

Upper mass limit for fit to GMC mass distribution 



a Restricted to logM„„-„. r / 6 [3.75, 3.8, 3.85, 3.9] for fits to the cluster mass distribution in the ring region. 

The slopes obtained for the differential mass distributions constructed using different binning strategies (i.e. the black and filled 
grey histograms in Figures 17 and 18 1 are generally in good agreement for the range of n \,j and n\, m values that we consider, 
but tend to show a larger dispersiorTtnan the fits obtained from the cumulative mass distributions. For some environments, the 
distribution of slopes is especially broad, or shows evidence for bimodality. This is most evident for environments where there 
is a relatively sharp bend in the mass distribution within the mass range that we use to estimate the slopes, e.g. GMCs in the bar 
region (panel [a] of Figure 17 1, and clusters in t he m aterial arm (A2, panel [f] of Figure 18 1. Though we simply use the median 
of each environment for our analysis in Section 5.3 and capture the broader distribution in the uncertainties, it is worth noting 
that the relevant value of the slope in these environments may depend on the physical process under investigation, and range of 
masses over which that process is likely to be acting. 



30 



, r, H I I | I I I I | I I I I | I M M I I | I I I I | I I I I | I h 

[b] ring 



GO [a] bar 

l3 

•j— i 
$H 0.5 

H 



° 0.0 
1.0 

O 

• i— i 
^— » 

o 

Ctf 0.5 

u 

0.0 



M 



'I'll''' 




Ji ii i I I 



-3 -2 -1 

-J_ I I | I I I I | I I I I | I h 

[c] Alin 



''I''' 1 '! 



I I I M-.I I I 



-3 -2 -1 

-J_ I I | I I I I | I I I I | I M i n 

[d]Alout : 



' ' ' ' <■' 



w_ 



0.5 



.A,- 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

[e]Al 



Jl 



'■ 



i_L 



'T+" ' 



J_ I I | I I I I | I I I I | 

- 1 m A2 



-1 



'''I i-T T l i" 1 l.j_i.L 

-3 -2 -1 



I I I I I I I I I I I I I I l_ 
[g] up 



i — 



■' ' iLl i i i i I i i 



0.0 

I I I I I I I I I I I I I H 1 r, 

[h] down 

-0.5 



Lj_l. 



^" I ' ^ 0.0 



-1 



Slope of GMC Mass Spectrum 



Fig. 17. — Histograms showing the distribution of slopes that we obtain from the trial fits to the cumulative (dot-dot-dashed lines) and 
differential mass distributions of GMCs in different M51 environments. The black histograms indicate the slopes obtained from distributions 
constructed using bins of equal logarithmic width. The filled grey histograms represent the distribution of slopes for a mass function constructed 
using an equal number of GMCs per bin. 
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Fig. 18. — The same as FigurefTT] but for young stellar clusters. 



